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SUMMARY 


This report is a user’s manual for version 2.1 of pEst, a FORTRAN 77 computer program for interactive 
parameter estimation in nonlinear dynamic systems. The pEst program allows the user complete 
generality in defining the nonlinear equations of motion used in the analysis. The equations of motion 
arc specified by a set of FORTRAN subroutines; a set of routines for a general aircraft model is 
supplied with the program and is described in the report. The report also briefly discusses the scope 
of the parameter estimation problem the program addresses. The report gives detailed explanations 
of the purpose and usage of all available program commands and a description of the computational 
algorithms used in the program. 


INTRODUCTION 

Parameter estimation techniques, in one form or another, have been in use at NASA Ames Research 
Center’s Dryden Flight Research Facility (Ames-Dryden) and other research organizations for many 
years. High-speed digital computers were first used for parameter estimation in 1968 (ref. 1), and 
the number of parameter estimation computer programs available has since greatly increased. The 
MMLE3 program (modified maximum likelihood estimation program, version 3) developed at Ames- 
Dryden (ref. 2) has been accepted as an industry standard for aircraft parameter estimation and has 
been used on a variety of aircraft programs. The MMLE3 program is representative in two respects 
of the majority of parameter estimation programs currently in use. First, it is designed for use solely 
in a batch processing environment. Second, the equations of motion defining the dynamic model used 
in the program are linear. For a large class of well-behaved parameter estimation problems, these two 
characteristics pose no serious limitations in the utility of the program. 

Recent flight test experience at Ames-Dryden has pointed out some of the limitations inherent 
in current parameter estimation programs. The dynamic behavior of aircraft at the extreme flight 
conditions currently being explored often cannot be appropriately modeled using the simple linear 
dynamic equations of motion. More accurate and flexible nonlinear models are often needed. The 
difficulties associated with extreme flight conditions, as well as those associated with the unique aircraft 
configurations currently being flown, have required significantly more attention from the analyst than 
previously. Interaction between the analyst and the estimation program is often the only viable means 
for obtaining results in a finite amount of time. 

In response to these problems, Ames-Dryden researchers have developed anew parameter estimation 
program, pEst. The pEst program is designed to be fully interactive; however, it can be run in a batch 
mode. The program supports full nonlinear capability in the dynamic equations of motion; linear 
equations are acceptable as a subset. 

This report documents the design philosophy, capabilities, and operational use of the pEst program. 
Section 1 defines the parameter estimation problem that pEst solves. Section 2 describes the philosophy 
of interactive program design as implemented in pEst. Section 3 describes the external files used by 
the program and how the program interfaces with other p 'grams. Section 4 gives a description of each 
command in the pEst command set. Section 5 discusses various algorithms available during program 
use. Section 0 defines the standard user routines supplied with the program. Both the equations of 
motion and the definition of all system variables used in the equations are included. In this manual, 
file names, program prompts, and literal program input are shown in italics to distinguish them from 
other text. The appendixes contain information on file formats used by the program (app. A) and 
listings of the help files used in the program (app. B). 


1 THE PARAMETER ESTIMATION PROBLEM 

Conceptually, the parameter estimation problem is straightforward: We are studying a physical system 

°< of motion that (hopefully) dsLL a mod J of t£ 

actud system^ We presume to know the form of the equations but not the values of certain parametric 

.rtu^.Tot.T^e v* P * rf ° m “ *“» U. ottual physical 

the input to the tytem and the reeponae of the system to the input. We seek to Infer the values of 

m«Ju”ed 0 ZponH“‘ t "‘ * "“' ir ^ “ ' h ' m0 ‘ W “ tU *«'•« -i‘k «h« 

_.. The . pE ^‘ P '°‘ r * m d °“ ,wo F “«. ptopam defines quantitatively the criterion for 

measuring the apoement between the model's computed response and the measured response. Second 
it mechanizes the search for the unknown parameter values. ’ 

, F ‘ gUre l lUustrates th * PEst parameter estimation process. The number in each block refers to the 
section in this manual describing the function of the block. 
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Uncertainty bounds 


Figure 1. The pEst parameter estimation process. 


1.1 Cost Funct : .on 

J''? "! le ' i0 " i * 1 •“■"tost function that is an explicit function of the computed response and thui 
an implicit function of the vector of unknown parameters. The cost function used in the program is 

<1 fit 


1 


/,{) = 2^7 E M»i> - i(Mrw'Wti) - iff,)] 
M 1 1*1 


whore n, and n, are the numbers of time history points and response variables respectively t is the 
ime variable, W the response weighting matrix, x the measured response, x the response computed 
y m egrating the equations of motion, ( the parameter vector, and superscript * denotes transpose. 
The cost function is quadratic in the computed response i. 
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1.2 Equations of Motion 

««-* differentia] ^nations 

response equation: ? th€ continuou8 -““e »tate equation and the discrete-time 
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where x is the corrected estimate of the state variable and X is the predicted estimate 

pni^. nZH ™“ZZs'Zl" ?r y The tCSP °" SCS ~ d “ •»* sample 

.he response function /Znole^unZn, T'T ^ ^ / Z* 

form of the functions / and a is defined bv a spt f ^ i Statei * nput ’ time- The specific 

I’ 1 '"' 1 witl > 'he program are documented in sectioTs'l” 10 ' * * SUbr0utinet ir ‘ P^'i the equations 
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All program commands and capabilities are aShtefor 7^7' an t imV’ a " d ' d r ' v “ 11 ‘""'face 
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simple and responsive mterfaceZZrZvM thfuar Z I J r ° eram ' ‘".''active or otherwise. A 
mistakes, In the interactive environment’ . n j . j,. 6a y ,ncr eased opportunities for making 

h«v« implemented error detection at everypotenti^rir 4 ^ increasin «'y im P ort ^t- We 
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program. Each program command has its own help file detailing the use and syntax of the command. 
Additionally, we have incorporated information on subjects of interest to the program user into the help 
file system. All help files are accessible during program use; they are also available outside the program. 

Interactive programs are generally reserved for difficult problems, where the approach to the solution 
is not clear at the outset. The solution progress tends to be discontinuous and incremental, with 
numerous dead ends met during the process. Efficient interactive problem solution requires a means 
of recovery from such dead ends. We have implemented a program feature that greatly enhances the 
potential to recover from errors and to restart the program if necessary. We have defined a program 
status file on which you can record the current status of the program at any point, thus allowing you 
to maintain an ancestor that can be used in the event of reaching a dead end. Effective use of the 
status file gives you freedom to experiment with different approaches to the problem, without fear of 
losing any progress already made. 

The pEst program is an interactive program. Some problems do not require much user interaction 
to obtain a solution. If you can define a sequence of pEst commands that, when executed, will solve 
the problem, you can use the program in a batch operating mode. 


3 INTERFACE TO OTHER PROGRAMS 


Tin* pEst package must be installed on your computer. Depending on the operating system and specific 
installation, a few features of the program may not be available (notably the help facility). 

The pEst package consists of three separate programs. The pEst program itself is the parameter 
estimation program. The thPlot program plots time history data and time history fits. The GetData 
program (ref. 3) selects signals and maneuver times for analysis. In use, the pEst program interfaces 
with these and other programs through several files external to the program. All file names used 
by pEst and other programs are italicized in the following sections only to distinguish them from 
other text. 

3.1 Measured Time History Data 

A file of measured time history data for the case to be analyzed must be available in a form suitable to 
the program. The program reads the entire measured time history file into memory; the upper limit on 
the number of time points that the program can handle is 2000. This limit can be easily changed by 
a single-line modification in the program code. While you can interactively select subsets of the time 
history for use in the analysis (see section 4.7.8), it is most convenient for the measured time history 
to approximate the time interval or intervals to be used. The time between sample time points on the 
file need not be constant. The measured time history data file is never rewritten or altered by pEst. 

Operationally, the measured time history file is divided into one or more time intervals called ma- 
neuvers. Each maneuver is treated as a complete and separate time history record in the integration of 
the equations of motion; the integration is reinitialized at the beginning of each maneuver. All maneu- 
vers are used together in t b estimation process; program variables, including estimated parameters, 
apply to all time points in Maneuvers. The program automatically breaks the time history into 
maneuvers when reading the h a new maneuver is defined when one of two conditions is found. Both 
conditions arc based on time values for successive time points read. If the time is nonincreasing or if 
the time increment exceeds a certain value (1.0 sec), then a new maneuver is started. In the simplest 
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and most common case, the file consists of a single maneuver. The maneuver times defined by the 
program may be displayed (see section 4.7.8). J 06 

can h* f eXP€CtS the ! ime hi8tory data to be on a named measured. Other file names 

can be specified by user command. ames 

The measured time history file is normally produced by the GetData program which selects the 
desired times and signals from the available data. The pEst program recognL signk by^hefrtames 
i rtfore the signal names on the file must match those expected bv nFst ’ 

SfSS,"* Tr rr- ^ 

6 ' 2 ' The GetData pr °S ram Md ths Me format are documented separately in reference 3. 

3.2 Program Status File 

Ihe program uses a status file to store the operational status of the program With the eveentm r 
time history data, the status file stores the value of every program variable Ffferf ^ exception of 
Hie , eo„ tf a. to efficient use of ,h. program The 

n tz *¥• ■* “ 

conclusion of a run. y ' d ' ,l C3n St ° re sumn » a ry results at the 

Efficient production use of pEst requires a status file at program startup. This file is not strictlv 
<1 med, as all variables are initialized with default values, which you can then change interactively 
over, manually setting the many variables for each case is laborious if many cases are to be analyzed' 

The initial status file can be obtained from one of several sources. For a large project vou will 
normally want to write a program that automatically creates an initial status file for each cLr this 
program should get starting parameter estimates from the simulation data base A status ITffrn 
whatever source) for one case can be copied and used initialise 

siaZVrmmUatch. mteraCt,VOly ' ‘ h,! r ‘"’“ irC ' 1 m ° dm 'a«ons is likely to take less effort than 

Jpzzzz ts,- srj v : ziz 

=-: ~ ir-arsss s: ;sr ■ ~ 

N h'°" 1 [ Esl “ ,C! * file that re(l ' C,S thc con >P let<! P ro « ram statu. prior to tertnina. 

‘ ‘ ‘ f n be uscd h y an y Program that analyzes or displays the results of DEst Tt ran *l c h 

Wl11 ' »**•><>-* ntodif, cations a, an initial status file for later casl " 

Jllti. Pr ° gram “ PeC,S lhe W named ot, '« "“><» ran be specified by 

1 lie format of the status file is documented in appendix A. 

3.3 Computed Time History Data 

A fde Of liutr Ids, or, data contputed by the program is available for use outside ,h. pEst program 
also * h ' ‘h” 10 ' Pr ° g ' am fit' f, can 
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I lie format of the computed time history file is identical to that of the measured time history file. 
The default name for the file is computed ; other file names can be specified by user command. 

3.4 Command Files 


1 he pEst program is an interactive program, with commands normally entered singly from the terminal 
keyboard. Sometimes, however, you might have a sequence of commands that you will be executing 
as a group more than once. The program provides a way of automating such repetitive tasks. If you 
write the desired sequence of commands on a file, you can then instruct pEst to execute the commands 
from that file as though they were typed from the keyboard (see section 4.8.1). After executing all the 
commands in the file, control is returned to the terminal. 

”1 here is no default name for the command file; any name can be specified by user command. 

3.5 Plot Commands File 


The pEst program runs the thPlot program when plotting time history fits. A temporary scratch file 
is created to communicate information from pEst to thPlot whenever plotting. In normal operational 
use, this file is deleted after the successful completion of the plotting, so you should never be aware of 
its existence. In the event of catastrophic program failure, however, it is possible that this file, named 
phstJ.hPlot.temp , might remain in existence. 


4 HOW TO RUN THE PROGRAM 

The pEst user interface is command driven; there is no menu. There is a single level of interaction 
between you and the program; all commands are available at any time. The program will accept com- 
mands when it issues the prompt pest command. The program recognizes a large set of commands; it 
is your responsibility to know the available command set and enter an appropriate command. Each 
command starts with a simple English key word and is followed by optional specifications. An inter- 
active help facility is available to help you find both the appropriate command and the proper syntax 
of a command. 

Each command controls a different aspect of the program; the applicable set of specifiers and their 
syntax arc command dependent. Several specifiers, however, are used in the same syntax in several 
commands; their usage and syntax are documented in the following sections. 

Several commands use switches to specify program options. A switch has one of two values; 
it turns an option either on or off. Most switches have one or more synonyms and antonyms, which 
are documented in the help files. Which name you use is a matter of personal preference and context; 
some fit more naturally into a certain command context than others. An algebraic sign is part of the 

switch specification; inverting the sign inverts the value of the switch. For example, -false is equivalen* 
to +true. 

Several commands also use key word and value pairs to specify program options. A key word and 
value pair assigns a value to a program variable. The key word is a word or abbreviation recognized 
by the program that corresponds to a variable in the program. The key word is delimited by a space 
or an equals sign and is followed by a value of the correct type. 

All input to pEst is case insensitive; you can use any mixture of upper- and lower-case letters. 
Most key words, program variable names and values, and switches ran be abbreviated. The minimum 
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1 he following sections document the individual commends of the pEst commend set. 


4.1 Help Command 

comp^onsiTand del^^^ t 0 J7S^ in ^ P T ^ ^ ° nUne he,R pr0vides 

There are three classes of inform aH™ ?, prog comn >ands, variables, and subjects of interest. 

help commands, help voneWcI aLTel Z.c. - ^ *° pi “ The 

available commands, program' variables and tcfniTs^ 0 Stlngs an< | short (one-line) descriptions of all 
specific command, progrfm variable « l 77^7 V V ° U U8C ^ with the " a ™ of a 
on syntax and usage and a few examples wiieiT! ** \ *77 help file list,ng givin S information 
or variable name, you mu!t ule ThTfuU nr^a ^ T* / ^ With . a Spedf,C command 

abbreviations are not acceptable If you use the hel7 ° 1 6 ^° m ™ and or vanab le; synonyms and 
description of the help facibty P C ° mmand W ‘ th n ° arguments . you get a brief 

All help file listings are included in appendix B. 

ing ou the OKratk^em.TnTr^rw ^plemimabTe ™ sZlTy's-ems 0 " ^ ^ deP "" i ' 
Examples of the help command are 


HELP 

Help variables 
help iterate 
HELP PARAMETERS 


4.2 Program Startup Commands 

1 lie program expects two files to be available at start.,™. ™ , . 

initial status file. If they are available n F t n measured time history data file and an 

the file measured and the “tip™ rim tlZ fZ * ^ ^ tk 

tltc respective files Two command, are e e current Both are default file names for 

without exiting the prig, am *° g * 1 StartU|> dala f '“ m «'« or to restart 


gram' mctly. ~t n 7"Z pL^ T7" " ‘° * T“" * M~r daU file into 

sru lion 6.2) are found on the file If a signal' ^'001'^ "V *" ?' g "* ls dcr ‘ ,l ' (l in 1,10 P'ogram (see 
used for ,| 1C timc historv “ ”° l ,0 “ n<i rd «- » constant value of rero is 

ami windows (see Bcction 4 7 8) are automTt7H° 8ra 7 Var |! abIeS ’ < ’" Cfl “ n,aneuvers (see section 3.1} 
e«n specify tho name of romma " d " «*«-■ ^ 

;; ,iW ' " * ~ Thp P ro g rl| u automatically attempt 

cannot ■“* ™ d the «• *» „ 0 
command lino. err ° r me88age 18 P r *nted, and control is returned to the 
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If there are no measured time history data in program memory, some commands will not be accepted 
by the program. In particular, any command requiring integration of the equations of motion (such 
as iterating) will not be accepted; if you attempt any such command, an error message is printed, and 
control is returned to the command line. 

Examples of read commands are 

READ 

Read measured. case6 


4.2.2 Restore command. — Use the restore command to read the program status from a status 
file. You can specify the name of the file to be read; if you omit the file name, the program uses the 
name last specified in a restore command. The program automatically attempts to read the file current 
at startup. 

If the file does not exL t when you attempt to read it, an error message is printed, and control is 
returned to the command line. If the program fails when reading a file, all program variables successfilly 
read into program memory up to that point are retained, an error message is printed, and control is 
returned to the command line. 

Examples of restore commands are 

Restore 

RET»T current .f Itl8.man6 
4.3 Program Termination Commands 

Several commands are available for terminating the program. Program status can be saved at program 
termination if desired. 


4.3.1 Save command. — Use the save command to write the program status to a status file. 
You can specify the name of the file to be written; if you omit the file name, the program uses the 
name last specified in a restore command. The default file name at program startup is current. If a 
file with the specified name exists, it is overwritten. 

If the program fails when writing a status file, no status file is written, an error message is printed, 
and control is returned to the command line. 

Examples of save commands are 


Save 

save curr.linal 
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1.3.2 Quit command. — Use the quit command to write the current program status to a status 
file and then terminate the program. This is exactly equivalent to running, in order, the save and abort 
commands. Note that you cannot specify the name of the status file to be written; the program uses 
the name last specificed in a restore command. 

If the program fails when writing the status file, the file is not written, an error message is printed, 
and the program terminates. 




4 
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quit 


An example of the quit command is 


.« •“ U “ ‘iV 4 * “ n ““ <1 *•"*“* the Profiram and return control 

JZyZITXTZ ” ot “ ,e ,he curreBt piogram su, “ : ^ — 

An example of the abort command is 


Abort 


4.4 Plotting Commands 

rwlabta ' t ‘" PIOt P '°* r “ 1 fOT Pl0,ti ” S ,imC hiS, ° rieS of ^PO" 5 ' fits and 

toory d" u PES ‘ with thPlot through a hie containing computed time 

Th ® C | qUati °uf ° f m ° ti0n are inte 6 rated to Produce the time history of the computed variables All 


4.4.1 Write command.— Use the write command to write the computed time history data from 

:7:~n V°ZT tn ^ le - “ torie8 for •“ computed state 

t to the file. The time histories of all computed variables are made consistent with the current 
I gram status prior to writing the computed time history file by integrating the equations of motion 
“ 7 l " :tU ' re " , P ro 6 r am status. Time histories for all response residufl 

control, response, and extra variables, are also written to the Me. You can specify the name rf tte fib 

0 >c written; if you omit the file name, the program uses the name last specified in a write command 

1 ho default file name at program startup is computed. command. 

If a file with the specified name already exists, it is overwritten. If no measured time history data 
hxamplcs of the write command are 


write 

Write computed . maneuver lOg 


... r 4 :V 2 . Plot command -— Uso the £lot command to plot time history fits for all active response 
amble, (see section 4.6.4) using the thl-lot program. The computed time histories are first wbZ 
an external file by automatically running a write command. The time histories are plotted four per 
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page; multiple pages are plotted if required. The measured response and computed response are both 
plotted on the same axis; the axis is automatically scaled to accommodate both variables. You can 
optionally request plots of additional time history variables (response residuals, states, controls, and 
extras); the plots for these variables follow the plots of the time history fits of the response variables. 

If no measured time history data exist in program memory, an error message is printed, and no 
plots are made. If the integration terminates prematurely for any reason, an error message is printed, 
but plots are still made. If the program fails to write the computed history file for any reason, an error 
message is printed, and no plots are made. 

Examples of the plot command are 

Plot 

plot i-res ids ♦states 


4.4.3 thPlot command. — Use the thPlot command to run the thPlot program from within 
pEst. The computed time histories are first written to am external file by automatically executing the 
write command. The pEst program does not communicate with thPlot when using this command* you 
have complete freedom to read any files and plot any signals desired. After thPlot terminates, control 
is returned to pEst. 

An example of the thPlot command is 


thplot 


4.5 Iterate Command 

Us«- the iterate command to start the iterative parameter estimation process. You can specify the 
maximum number of iterations and the desired minimization algorithm. If you do not specify the 
number of iterations, no iterations are done; the equations of motion are integrated, and the cost 
function is evaluated. The three minimization algorithms available are gradient, Newton- Raphson, and 
Oavidon- Fletcher- Powell. If you omit the minimization algorithm, the algorithm last specified remains 
in effect. The default algorithm at program startup is Newton- Raphson. All the iterations specified 
with a single iterate command use the same minimization algorithm. 

Before attempting to iterate, the program tests the validity of the current program status for 
estimation; if the validity test fails, iterating is not allowed, and control is returned to the command 
line. For example, if no measured time history data are in program memory, iterating will not be 
allowed. If all validity tests succeed, the program iterates until either the specified number of iterations 
is completed or the convergence criterion (see section 4.7.4) is met. The iterative process may terminate 
wit h one or several error messages if any of numerous problems is detected. If the estimation is 
prematurely terminated, all progress made up to the point of failure is retained and is reflected in the 
current program status. 

Examples of the iterate command are 


6 newton-raphson 
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Itsrats 
it 3 
IT 


4.6 System Variable Commands 

The equations of motion define the dynamic system analyzed by the program. The equations of 
motion contain numerous time history and parametric variables defining the specific characteristics of 
the dynamic system. The system variable commands allow you to display and to modify these variables 
and therefore the characteristics of the dynamic system. The time history and parametric variables are 
vector variables, and each vector element has several characteristics or attributes. Each system variable 
command allows you to display and to modify the attributes of selected elements of a specified vector. 

The first part of each system variable command selects elements of the vector to be displayed or 
modified. The syntax of this specification is common to all system variable commands. A generalized 
list of vector elements follows the command name. The list can be a list of element names, delimited 
by commas or blanks. It can also be one of two key words, qll or active ; all selects all elements of 
the vector, and active selects the vector elements with active status. The definition of active status is 
dependent on the vector and is defined in each of the following command descriptions. If the list is 
omitted, all active elements are selected by default. Whenever you execute a system variable command, 
the values of the selected elements are displayed. 

For each system variable command, several optional descriptors control the attributes to be 
displayed or modified. The order in which they are specified on the command line is immaterial. The 
descriptors and the syntax for using them are documented in each of the following 
command descriptions. 


4.6.1 Parameter command. Use the par ameter command to display and modify the at- 
tributes of the parameters. The parameters are the variables that can be estimated by the program. 
Each parameter has five attributes: current value, estimation status, predicted value, Cramer-Rao 
bound, and change in value from previous iteration. The current values of the selected parameters are 
always displayed after the command is entered. The estimation status of a parameter is active or inac- 
tive, ,m active parameter is one that is currently being estimated, while an inactive parameter is one 
that is not. You can modify the estimation status of the selected parameters with the ±active switch. 
1 lie predicted value of a parameter is fixed and can only be displayed with the + predicted switch. The 
Cramer-Rao bound and the change in parameter value from the previous iteration are defined in the 
estimation process and cannot be directly modified by the user; they can only be displayed with the 
±bound and ±delta switches, respectively. You can modify the current value of a parameter in one 
of two ways. You can explicitly specify a value, and that value is assigned to all selected parameters. 
Alternatively, you can turn on the ± rest ore switch, and each selected parameter will be reset to its 
corresponding predicted value. 


The parameters defined by the standard user routines are documented in section 6.2.1. 
Examples of parameter commands are 


PARAMETER clp«-0.25 
par cma.cNorna cnda.cNorada *on 
Par Activa ♦Raatora 
parm *cr ♦delta 


4.6.2 Constant command. Use the constant command to display and modify the constants 
in the program. The constants are program variables that cannot be estimated by the program. The 
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value of each selected constant is always displayed after the command is entered. You can modify the 
value of a constant; the specified value is assigned to all selected constants. 

The constants defined by the standard user routines are documented in section 6.2.2. 

Examples of constant commands are 

constant aass«10S6 . 0 
COIST ixy.iyz 0. 

Const all 


4.6.3 State command. — Use the state command to display and modify the attributes of the 
state variables in the program. A state variable has two attributes: its status and its integration limit. If 
a state is active, the state equation for the state is integrated, and this integrated value is used in the 
equations of motion. If a state is inactive, the state equation is not integrated, and the equations of 
motion are modified to remove the corresponding state equation. You can modify the status of the 
selected state variables with the 4- active switch. The integration limit for a state variable is used to 
avoid catastrophic program failure during integration of the equations of motion. If during integration 
the value of a state variable exceeds its limit, an error message is printed, and the integration is 
terminated. The integration limit is a floating-point value specified using the lim it key word and 
value pair. 

The state variables defined by the standard user routines are documented in section 6.2.3. 

Examples of state commands are 

state v , alpha, an, q +on 
STATE p r 1U-10000. 


4.8.4 Response command. — Use the respo nse command to display and modify the attributes 
of the response variables in the program. A response variable has two attributes: its status and its 
weighting. The status of a response variable determines whether or not the response variable time 
history is computed. If a response is active, the response variable time history is computed and made 
available for other uses. If a response is inactive, no response time history is computed. You can 
modify the status of the selected response variables with the + active switch. The weighting of a 
response variable specifies the variable’s weighting in the cost function and is specified with the weight 
key word and value pair. The key word outp ut is a synonym for response. 

Note that making a response variable inactive is not equivalent to making its weighting zero. If a 
response variable is made inactive, the equations of motion are modified to remove the corresponding 
response equation. There may also be secondary changes in the equations of motion to remove all 
usage of the response variable. These secondary changes depend on the equations of motion used. If 
the response is active but has zero weighting, the response is computed and used in the equations of 
motion hut does not directly influence the cost function. 

The response variables defined by the standard user routines are documented in section 6.2.5. 

Examples of response commands are 
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response beta p r ♦on 
OUT alpha *>160. 


4.8.5 Flag command. — Use the flag command to display and modify the flags in the program. 
A flag is a logical variable used in the equations of motion. The flags typically select alternative forms 
of the equations of motion or sources of data. You can modify the value of a flag with the -bon switch; 
the switch value is assigned to all selected flags. 

The flags defined by the standard user routines are documented in section 6.2.7. 

Examples of flag commands are 

flag use.avg.qbar 

Flag use.avg. alpha, uss.&vg.beta +0N 
FLAG all +off 

4.7 Program Variable Commands 

Use the set and show commands to display and modify the program variables and options controlling 
the estimation process. 

The set command sets the value of program variables. With one exception, the command syntax is 
t ho command key word followed by a key word and value pair specifying a program variable and setting 
its value. The valid key words and the program features they control are described in the following 
sod ions. Each key word is given with its first few letters underlined; the underlined portion is the 
minimum abbreviation recognized by the program. The value type is dependent on tbe variable; for 
each variable described, the range of legal values is specified. Any variable modified by the command 
is displayed. 

The show command allows you to display the values of program variables. The command syntax 
is the command key word followed by a variable name or appropriate abbreviation. In addition to the 
variables described by the set command, there is one variable that is only displayable. 


4.7.1 Integration method variable. — The integ Meth variable specifies the algorithm used 
when integrating the equations of motion. The variable has two possible values, euler and runge-kutta; 
the default value is runge-kutta. 

Examples of the use of the integration method variable are 

show integMsth 
SET intsg rungs 

4.7.2 Gradient method variable. — The gradM cth variable specifies the algorithm used when 
computing the finite-difference gradients. The variable has two possible values: single-sided specifies the 
forward difference algorithm, and Rouble-sided specifies the central difference algorithm. The default 
value is single-sided. 

Examples of the use of the gradient method variable are 

sh gradnsthod 
Sst GradM 2 
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4.7.3 Gradient delta variable. — The gradD elta variable specifies the parameter increment 
used in computing the finite-difference gradients. The parameter increment used for each parameter is 
defined by 

d( = gradDelta * max(£, 0.000001) 
where £ is the parameter value and is the parameter increment. 

Valid values are floating-point numbers; the program default value is 0.0000001. 

Examples of the use of the gradient delta variable are 

sh graddelta 
SET GRADD»0.0001 


4.7.4 Convergence bound variable. — The bound variable defines the convergence criterion for 
the estimation process. If the percentage change in cost between two successive iterations drops below 
the value of the bound variable, convergence is declared, and the iterative process is terminated. Valid 
values for this variable are floating-point numbers; the program default value is 0.0001. 

Examples of the use of the convergence bound variable are 

SHOW B0UXD 
set bound 1.0e-6 


4.7.5 Message level variable. — The ms g Level variable controls the amount of output that the 
program prints during use. Higher values produce larger amounts of output. Valid values for the 
variable are integers between 0 and 100; the default value is 50. The help file lists the significance of 
the various numerical values. 

Examples of the use of the message level variable are 

sho nag 

Set MsgLev-65 


4.7.0 Plot title variable. — The title variable specifies the title used on the time history plots 
of response fits. Valid values are character strings of up to 40 characters. If there are blanks embedded 
in the string, you must put the whole string inside quotes. The default value for the title is blank. 

Examples of the use of the plot title variable are 
Show Tit 

SET title 'Space Shuttle Flight 4 Maneuver 3b' 


4.7.7 Statistics variable. — The stat istics variable contains sample statistics of all measured 
variables from the time history data file. The statistics are defined whenever a measured time history 
file is read and may not be modified during use; they may be displayed only with the show command. 

An example of the use of the statistics variable is 

Sho Stats 
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The Jffimjou, variable specifies the window or windows used in the analysis. The window variable 
specification has three elements: the window number, the maneuver number, and the time ,^8^ 
Whenever a window is referenced, all currently dehned maneuvers and window, 1 dtaplly^ 
Examples of the use of the window variable are 


Show windows 
SET WINDOW TIME 
set wind 2 mam 2 


0-10 

time 0.5 7.5 


4.8 Advanced Commands 




Examples of the do command are 


do initialize 
Do Startup.X29 


Examples of system commands are 


SYSTEM FILES 
ays help copy 

ays Rename Current. init.f 18 Current 


5 ALGORITHMS 


The 

The 




5.1 Minimization Algorithms 


I ho pEst program has three algorithms available 
(or Rtcepcst descent), modified Newton- Raphson, 


for iteratively minimizing the cost function: gradient 
and Davidon-Fletclier- Powell. The Newton-Raphson 
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5.2 Gradient Computation 
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parameter increment used in computing the gradients is a fixed percentage of the parameter value (see 
section 4.7.3). 

Dotli single-sided and central difference algorithms are supported. Gradient computation using 
the single-sided difference algorithm requires n p + 1 solutions of the equations of motion, where n p is 
the number of parameters; the central difference algorithm requires 2n p solutions of the equations, of 
motion. The gradient computation requires a large number of solutions of the equations of motion; 
the computational time involved is typically the most significant percentage of the total time used in 
an iteration. The single-sided difference algorithm is the faster of the two algorithms, but it is also 
the less accurate. The increased accuracy of the central difference algorithm, however, is not really 
consequential until close to convergence; we generally find it most efficient to use single-sided differences 
for all but the last few iterations. The Gauss-Newton approximation (ref. 2) is used to compute the 
second giadient, thus the computational burden of computing the second gradient in addition to the 
first gradient is not significant. 


5.3 Integrating the Equations of Motion 

There are two algorithms available for integrating the equations of motion: forward Euler and fourth- 
order Runge-Kutta integration. The forward Euler integration requires one evaluation of the state 
derivative function / for each time point of the solution, while the Runge-Kutta algorithm requires 
four such function evaluations per time point. Both algorithms require a single evaluation of the 
response function g for each time point of the solution. Thus, the Euler algorithm is the faster of the 
two. The Runge-Kutta algorithm, though slower, is more accurate and has a larger region of stability. 

The nonlinear character of the functions / and g requires consideration of several issues that are 
not relevant for linear equations of motion. The inherent possibility of singularities in / and o means 
that catastrophic termination of the program (with subsequent loss of anything done up to that point) 
is a real consideration. Extreme caution (maybe even paranoia, as some have suggested) is necessary to 
anticipate and eliminate these possibilities. Singularities, however, are just extreme and obvious cases 
of ill conditioning; the more subtle cases can also generate equally catastrophic errors. In the estimation 
process, it is not uncommon for an intermediate solution to diverge greatly from the final converged 
solution. The large magnitudes of the state and observation variables in these intermediate solutions 
can easily exceed the bounds of validity of intrinsic FORTRAN functions. We have implemented limit 
chocking on the state variables (see section 4.6.3) to preclude catastrophic computational errors. 

6 STANDARD USER ROUTINES 

1 lie pEst program is supplied with a set of equations of motion to model a wide range of aircraft prob- 
lems. The equations of motion are a full six- degree-of- freedom nonlinear set of differential equations. 
We do not assume that the aircraft is symmetric. We do assume fixed aircraft geometry and constant 
mass characteristics. No propulsion or rotating-mass effects are included. We assume a flat earth and 
constant gravitational acceleration. We use English units throughout the equations. 


0. 1 Equations of Motion 

The standard user routines define the equations of motion used by the program. The equations are 
divided into the state equations and the response equations. All time history and parameteric variables 
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used in the equations are described in section 6.2. There is no explicit division of the equations 
into longitudinal and lateral-directional subsets. The state equation for airspeed is not implemented 
in pKst. 

The standard user routines define 8 state equations, 5 feedback equations, and 14 response equa- 
tions. Only in rare instances, however, will you use the complete set of equations. In practice you will 
probably interactively define a subset that is dependent on the maneuver being analyzed. 

You can independently activate the state equations and response equations. Each state equation 
is integrated only if it is active (see section 4.6.3); inactive state equations are not integrated and are 
removed from the equations of motion. Each response equation is evaluated only if it is active (see 
section 4.6.4); inactive response equations are removed from the equations. 

Each state variable on the right-hand side of the state and response equations can come from one of 
three sources. First, the computed time history of the state variable can be used, if available. Second, 
the corresponding measured response variable can be used, if available. Third, for some state variables, 
a constant value can be used. The program flags and the status of the state variables determine the 
source of each state variable on the right-hand side of the equations. Five flags specify using average 
values of measurement variables ( V , a, /?, 0, and <j>) in the equations. If a flag is turned on, the value 
of t he corresponding constant variable is used in the equations. If the flag is not turned on (or if there 
is no corresponding flag for the state variable), the status of the state variable determines the source. 
If the state variable is active, the computed state time history is used. If the state variable is inactive, 
the corresponding measurement variable is used. 

Two additional flags specify the source of the extra time history variables (such as q) on the right- 
hand side of the equations. If the flag for an extra variable is turned on, the value of the corresponding 
constant variable is used; otherwise the measured time history of the extra variable is used. 


6.1.1 State equations. — The predicted state variables x are obtained by integrating the state 
equations defined in the following equations. 

qsR 


where 


a 


= q - tan /?(p coso + r sin o) — 


] T p — Ixy9 Ixt* 
I y<) ~ IyiT ~ IxyV 
IzT ~ fxzP~ Iyzi 
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mV cos/3 


Cl 


+ 


g R 


V cos/? 


(cos 0 cos <f>cosa + sin 9 sin a) 


4> = 


psinor — rcosa + — — Gy 
mV 

4- ^[cos /? cos 6 sin <j> - sin /J(cos 6 cos 4> sin a — sin 9 cos a)] 

qsbCtR + [qr(Iy - /,) + {q 2 - r 2 )I vt + pql xt - pr/ ry ]/R 
qscC n R + \pr(I t - I x ) + (r 2 - P 7 Vxz + qrhy ~ P?f v *)/ R 
qsbC n K + \pq(I x - I y ) + (p 2 -- q 2 )I IV + prl yz - qrl xx ]/ R 
q cos <(> - r sin <t> 
p + tan 9(r cos 4> -I- q sin <f>) 
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is reference span, 
reference chord, 
coefficient of lift, 
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coefficient of rolling moment, 
coefficient of pitching moment, 
coefficient of yawing moment, 
coefficient of lateral force, 
gravitational acceleration, 
are moments of inertia, 
cross products of inertia, 
is mass, 
roll rate, 
pitch rate, 
dynamic pressure, 
yaw rate, 

conversion factor (57.2958), 
reference area, 
total velocity, 
angle of attack, 
angle of sideslip, 
pitch attitude, and 
roll attitude. 

6.1.2 Initial conditions.— The initial conditions of the state variables used in integrating the 
equations are defined as follows. 


V(0) 

= 

V io -V h + Vo 

a(0) 


(a*o - ai)/k a + 00 

m 

= 

{Pzo ~ Pb)/ k 0 + 0o 

p(o) 

= 
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/(« is the value of the state variable at the beginning of the integration, subscript b denotes 
„ * bias for the corresponding observation variable, subscript 0 denotes the increment to 

11: S state value, and subscript r„ denote, the measured value of the corresponding observation 

variable of the initial point. 


S.1.S Total force and moment coefficient.,- The total force and moment coefficients used 
in the state and response equations are defined as follows. 

Ci = C/vcosa - Casino; 

C A = C Ao + C Aa a+ j~CA'q + CA t '6' + CA' l 6i + CA h h + CA li k 
Cn = C No + C Na a+ -CnJ + C Ntt *. + C N , X 6 t + C Sfj h + C Nh h 
Cy = Cy 0 + C Yff 0 + £^(CY„P + Cy r r) + Cy u 6* + Cy tj 6 r + Cy h 63 + C Y(i 64 
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where 
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is coefficient of axial force, 
coefficient of normal force, 
aileron deflection, 
elevator deflection, 
rudder deflection, 
control surface 1 deflection, 
control surface 2 deflection, 
control surface 3 deflection, 
control surface 4 deflection, 

are axial force parameters, 

normal force parameters, 

lateral force parameters, 

rolling moment parameters, 

pitching moment parameters, and 

yawing moment parameters. 


All parameters are defined in section 6.2. 


fl.1.4 Response equations. — The computed response variables z are obtained by evaluating 
the response equations defined in the following equations. 
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is roll acceleration, 
pitch acceleration, 
yaw acceleration, 

are axial accelerometer position parameters, 
lateral accelerometer position parameters, 
normal accelerometer position parameters, 
aircraft center-of-gravity position constants, 
angle-of-attack measurement parameters, 
angle-of-sideslip measurement parameters, 


and the subscript z denotes the computed response value. 


6.1.5 State feedback equations. — The corrected state variables x are obtained by evaluat- 
ing the discrete feedback equations (see section 1.2) defined by the following standard user routines. 
Feedback is implemented for only five of the state variables. 

a = a + (g a /k a )(ai - a t ) 

P = 0 + (9filh)(fo-P.) 

P = P + 9p(Pz-Pz) 

7 = q + 9q(<lz - Qt) 
r = r + g T (r ; -r z ) 

where g n ,!l0,9p, 9q, and g r are feedback gain parameters, the superscript ' denotes the predicted state 
value, and the superscript * denotes the corrected state value. 


6.2 System Variables and Names 

'Fho standard user routines define the names of all system variables: parameters, constants, states, 
controls, responses, extras, and flags. In program use, all system variables are accessed by their 
given names. 


6.2.1 Parameters. — The standard user routines define 97 parameters. Parameter subsets are 
d<Tnicd below. Each parameter is individually documented in the parameters help file. 

6.2.1. 1 Stability and control derivatives: Any parameter whose name starts with the letter c is 
a stability or control derivative. There are three classes of stability and control derivatives: aerody- 
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namic biases, control derivatives, and stability derivatives. Each derivative is nondimensionalized; the 
details of the nondimensionalization depend on the parameter. All aerodynamic biases are completely 
normalized; they are dimensionless. All rotational rate derivatives are expressed in reciprocal radians. 
All o, f), and control derivatives are expressed in reciprocal degrees. The reason for the apparent 
inconsistency in nondimensionalization is historical and follows existing conventions. 

The stability and control derivatives can be divided into longitudinal and lateral-directional deriva- 
tives, though there is some overlap in the control derivatives. Tables 1 and 2 tabulate the longitudinal 
and the lateral-directional stability and control derivatives, respectively. (In the tables and text, math- 
ematical variables in parentheses follow the corresponding program variable names.) 


TABLE 1 — LONGITUDINAL STABILITY 
AND CONTROL DERIVATIVES 


Sta te or control 
Aerodynamic bias 
Angle of attack 
Pitch rate 
Elevator deflection 
dl deflection 
d2 deflection 
d3 deflection 


Axial 

force 

caO (C Ao ) 
caa (C Aa ) 
caq (C\) 
cade (C A(t ) 
cadi (C A(l ) 
cad 2 (C A J 
cad3 {C A ) 


Force or moment 
Normal 

force 

cNormO (C/v 0 ) 
cNorma (C^ a ) 
cNormq (CVJ 
cNormde (Cj^ f ) 
cNormdl (Cn 6 ) 
cNormd2 ( Cjy * ) 
cNormd3 {C^ t ) 


Pitching 
moment 
cmO (C mo ) 
cma (C m J 
CI ^q {Cm q ) 
cmde (C m «J 
cmdl (Cm h ) 
cmd2 (Cm 67 ) 
cmd3 {Cm ) 


TABLE 2 — LATERAL-DIRECTIONAJ STABILITY 
AND CONTROL DERIVATIA ES 


Force or moment 


State or control 
Aerodynamic bias 
Angle of sideslip 
Roll rate 
Yaw rate 
Aileron deflection 
Rudder deflection 
d3 deflection 
d4 deflection 


Lateral 
force 
cyO (Cy 0 ) 
cyb (C Yp ) 
cyp ( c Yp ) 
v'yr (C Yr ) 
cyda (C\ 
cydr ( C Y(t 
cyd3 (C Yi3 
cydd (C Y , 


Rolling 
moment 
“clO (C fo ) 
clb (C tt ) 
dp (C<„) 
dr (C tr ) 


Yawing 
moment 
TnO (C„ 0 ) 
cnb (C ng ) 
enp ( C np ) 
enr (C nr ) 


,) dda (QJ enda (C ng J 
) cldr (C e J endr (C n J 
) cld3(0 <3 ) cnd3 (C n<3 ) 

jcldUQ ) cnd4 (C n , ) 


(>■^.1.2 State initial conditions: Each state variable has a corresponding state initial condition 
increment, parameter. The initial condition increment is added to each measured initial condition The 
dimensions of each initial condition parameter are the same as those of its corresponding state variable 
State initial condition parameters are vO (V 0 ), alphaO (o 0 ), qO ( q 0 ), thetaO ( 0 o ), betaO (/3 0 ), pO (p 0 ), rO 
(**■»), and phiO (<t > 0 ). v ^ u/1 
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6 2 1.3 Instrumentation parameters: Numerous parameters characterize instrumentation installa- 
tion a.ul calibration. Each response variable has a corresponding response b.w^ 
value is added to the computed response variable value at each time point. The dimensions of 
response bias parameter are the same as the dimensions of the corresponding response variable. Re- 
sponse bias parameters are vBias (V b ), alphaBias (ob), qBias (ft), thetaBias (ft,), anBias KJ , ***** 
(« r j. qdotBias (ft), betaBias (A>), pBias (pt>), rBias (a), phiBias (<h>), ayBias (o^), pdotBias (pt), 

an<l rdotBias (rb). 

Several response-measuring instruments have parameters describing their location in the *rcraft. 
Each instrument-position parameter specifies the location of a response- measuring instrument aft of 
to the right of, or above the aircraft reference point. Instrument- position parameters correct computed 
response values for instruments not located at the aircraft center of gravity. In operation, the program 
determines the distance from the instrument to the center of gravity by subtracting the instrument 
position from the center-of-gravity position (see section 6.2.2.1). The dimensions of all instrument- 
position parameters are feet. Table 3 gives the instrument-position parameter names. 

TABLE 3 - INSTRUMENT POSITION PARAMETERS 


— 


Aircraft axis 


Response variable 

X 

Y 

z 

Velocity 
Angle of attack 
Angle of sideslip 
Axial acceleration 

xv (xv) 
xa (x a ) 
xb (x 0 ) 
xax (x 0 J 

yv (yv) 
ya (y a) 
yt> (yp) 

yax (y a> ) 

zv (xv) 
za (z Q ) 
zb (z/j) 
zax (z 0l ) 
zay ( 2 a> ) 
zan (z a .) 

Lateral acceleration 

xay (x„,) 

yay (»«,) 

Normal acceleration 

xan (x a J 

yan (y«.) 


Two parameters define the flow amplification factors for the aircraft flow angle sensors. The upwash 
factor for the angle-of- attack sensor is ka (fc a ), and the sidewash factor for the angle-of-sideslip sensor 
is kb (bfl). Both parameters are dimensionless. 


6 g, I Feedback gains: Five parameters define feedback gains for five of the state equations Valid 
values for feedback gains are between 0 and 1; a value of d implies no feedback^ While the feedback gains 
are parameters and can therefore be estimated, such use is experimental. Feedback gain parameters 

are gAlpha ( g n ), gBeta (y*?), gP (y p ), gQ (»«)> and (3r)- 


6.2.2 Constants. — The standard user routines define several constants. The use and dimensions 
of each constant are defined as follows. 


r , > •> i Aircraft physical characteristics: Constants defining the mass characteristics of the aircraft 
are miss'd"), ix (/,), iy </ v ), i* (/.), bey (/„). ™ (/«). and iyz (/„,)• The mass is * 

inertias in slug-feet squared. Constants defining the reference dimensions of the aircraft are area (s), 
span (fc), and chord (c); the area is in feet squared and the span and chord are in feet. Constants 
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defining the location of the aircraft center of gravity are xcg ( x e% ), ycg (y cf ), and zcg (z cg ); they are 
measured in feet aft of, to the right of, and above the aircraft reference point, respectively. 

6. 2. 2. 2 Time history variable averages: Several constants specify average values of measured time 
history variables. The avg_v, avgjdpha, avgJ>eta, avg-lheta, and avg_phi constants contain average 
values for the corresponding response variables. The avgjqbar constant contains the average value of 
tho dynamic pressure. Whenever a measured time history file is read (see section 4.2.1), the average 
value of each time history variable is computed, and values are assigned to the corresponding constants. 
Of course, at any point in the program, you can modify the value of a constant (see section 4.6.2). 
The equations of motion use the value of the constant variable containing the time history average of 
a response variable if its corresponding flag is turned on (see section 6.2.7). 

6.2.3 States. — The standard user routines define eight state variables. All state variables are 
used "internally in the program in English units. The wind-relative velocity is defined in spherical 
coordinates by three states: airspeed, angle of attack, and angle of sideslip. Airspeed is measured in 
feet per second, and its state name is v (V). Angles of attack and sideslip are both measured m degrees 
and are named alpha (a) and beta (/3), respectively. The aircraft rotational velocities are defined m 
the aircraft body axes and are all measured in degrees per second. The roll, pitch, and yaw rate states 
are named p (p), q (9), and r (r), respectively. The aircraft attitude is defined by the the aircraft Euler 
angles, measured in degrees. The pitch and bank angles are named theta (6) and phi ( <t > ), respectively. 
Heading angle is not used in the aircraft equations of motion. 


6.2.4 Controls. — The standard user routines define seven control variables. All control variables 
are measured in degrees. The three conventional aircraft control surface deflections, elevator, aileron, 
and rudder deflection, are named de (tf e ), da {6J, and dr (6 T ). The program defines four additional 
controls, dl (61), d2 (i?), d3 (63), and d4 (64). 


0.2.5 Responses. — The standard user routines define 14 response variables. Each of the eight 
state variables (see section 6.2.3) has a corresponding response variable; the name and dimensions of 
each response variable are identical. Six additional response variables have no corresponding states. 
The aircraft rotational accelerations are defined in the aircraft body axes and are measured in degrees 
per second squared. The roll, pitch, and yaw accelerations are named pdot (p), qdot ( q)> and rdot (r), 
respectively. The aircraft linear accelerations are defined in the aircraft body axes and are measured 
in g. The axial, lateral, and longitudinal accelerations are named ax (a r ), ay (a y ), and an (o n ), 
respectively. 


0.2.6 Extras. — The standard user routines define three extra signals: Mach number, dynamic 
pressure, an d altitude, named mach (M), qbar ( q ), and alt (A), respectively. 

6.2.7 Flags. — The standard user routines define seven flag variables. The flags select alternative 

forms of the equations or sources of data. Each flag specifies using the average value of a measured 
time history variable in the equations of motions in place of a time- varying quantity. This option is 
particularly useful when you do not have a measurement time history for a particular variable. In 
surh a case, you can input an average value using the appropriate constant (see section 6.2.2) and then 
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the flag to u« tie average «1«. ia «h. equation.. The lag. are named ua.-a.g-., uae.avg.beta, 
uscJtvgJtheta, use ^vg 4) hi, use-avgjnach, and usejwg-qbar. 


National Aeronautics and Space Administration 
Ames Research Center 
Dryden Flight Research Facility 
Edwards, California , September 9, 1986 
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APPENDIX A— PROGRAM STATUS FILE FORMAT 


The program status file is a FORTRAN formatted file that stores the status of every program variable. 
Each record on the status file corresponds to an individual program variable or option; the records are 
grouped together into common record types. When the file is written by the program, a complete file 
is written; all records defined for the file are written. If you create the file outside the pEst program, 
it is not necessary to fill in all fields on the file. At startup, the program assigns a default value to 
each program variable; the value is retained unless it is redefined by reading a corresponding record on 
the status file. Except for the first record, which specifies the file format, the order of the records is 
immaterial. Any record can be repeated; if repeated, the last record read overrides all previous records. 
All fields are case insensitive; upper and lower case can be mixed. 

Each record type is defined by its key word; the key word is the first field of each record and is 
loft -justified in the first eight columns of the record. Each key word can be followed by several fields; 
the content and format of the fields, of course, depend on the record type. However, all fields of a 
common type have identical format, independent of the record type. There are no spaces between fields. 
Any field specifying a system or program variable name is 16 characters wide; the name is left-justified 
in the field. Any floating-point field is also 16 characters wide and is read with a FORTRAN gl6.10 
format. Any logical field is 2 characters wide, and contains either a t or an f, right-justified. All other 
fields are record dependent and are described in the individual sections that follow. All key words and 
literal phrases used in the status file are italicized in the following sections to distinguish them from 
other text. 


A.l Version Record 

The file has a single version record. The version record specifies the format used when reading the 
file, hence the version record must be the first record on the file. The key word for the version record 
is version. Four fields follow: The first is a character string with the phrase pest-current left-justified 
in a field of 16 characters. The program recognizes the file by this phrase, and it is required. The 
second field specifies the version number of the current file as a character string, left-justified in a field 
8 characters wide. The current value is 2.1. The third and fourth fields specify the date and time when 
the file is written by the program; these fields are not read by the program. 

An example of a version record is 

version pest-current 2.1 31-Mar-86 15:33:30 


A. 2 Title Record 

The file typically has a single title record. The key word lor the title record is title. A character field 
40 characters wide specifying the title follows the key word. 

An example of a title record is 

title Fi8 flight 11 long man 3 
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A.3 Parameter Record 


The file can have several parameter records; typically there is one for each parameter defined in the 
program. The key word for a parameter record is param. Five fields follow: The first is the parameter 
name. The second and third are floating-point fields specifying the predicted and current values of the 
parameter, respectively. The fourth is a logical field specifying the estimation status of the parameter. 
The last is a floating-point field specifying the Cramer-Rao bound for the parameter. 

Examples of parameter records are 


param cNormde 

param clp 


. 7798399770E-02 .8341009928E-02 T .0001045977 
-.3680106990 -.3680106990 F .0000000000 


A. 4 Constant Record 

The file can have several constant records; typically there is one for each constant defined in the 
program. The key word for a constant record is const. Two fields follow: The first is the constant 
name. The second is a floating-point field specifying the value of the constant. 

Examples of constant records are 

const avg. alpha 8.165277248 

const ixz 1870.199950 

const zcg 5.317500110 


A. 5 Flag Record 

The file can have several flag records; typically there is one for each flag defined in the program. The 
key word for a flag record is flag. Two fields follow: The first is the flag name. The second is the value 
of the flag, in a logical field. 

An example of a flag record is 

flag usa.avg.alpha F 


A.0 State Record 

The file ran have several state records; typically there is one for each state variable defined in the 
program. The key word for a state record is state. Three fields follow: The first is the state name. The 
second is a logical field specifying the status of the state equation. The third is a floating-point field 
specifying the integration limit for the state. 

Examples of state records are 

atata alpha t 100000.0000 

atata p F 10000.00000 
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A. 7 Response Record 

The file can have several response records; typically there is one for each response variable defined 
in the program. The key word for a response record is output. Four fields follow: The first is the 
response name. The second is a floating-point field specifying the average value of the measured 
response variable. The third is a logical field specifying the status of the response equation. The last 
is a floating-point field specifying the response weighting in the cost function. 

Examples of response records are 

output alpha 8.165277248 T 500.0000000 

output phi -.3456887860 f 10. 00000' 00 


A.8 Control Record 

The file can have several control records; typically there is one for each control variable defined in the 
program. The key word for a control record is input. Two fields follow: The first is the control name. 
The second is a floating-point field specifying the average value of the measured control variable. 

Examples of control records are 

input do -2.275935836 

input d3 .0000000000 


A.9 Extra Record 

The file can have several extra records; typically there is one for each extra variable defined in the 
program. The key word for a extra record is extra. Two fields follow: The first is the name of 
the extra signal. The second is a floating-point field specifying the average value of the measured 

extra variable. 

An example of an extra record is 

extra qbar 112.4635873 


A. 10 Maneuver and Window Records 

The file can have several maneuver or window records, or both; typically there is one foi each maneuver 
or window defined in the program. The formats of the two records are identical; the only difference 
is the key word. The key word for a maneuver record is maneuver , while that for window record 
is window. Two fields specifying the start and end times follow the key word. The format for each 
time is identical. There are two blank spaces, followed by a two-digit hours specifier, a period, a 
two-digit minutes specifier, a period, a two-digit seconds specifier, a period, and a three-digit milli- 
seconds specifier. 

Examples of maneuver and window records are 

maneuver 10.42.18.608 10.42.25.487 

vindov 10.42.18.608 10.42.25.487 
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A. 11 Option Records 

'Die file can have several option records; typically there is one for each program option defined. The 
key word for an option record is option. Two fields follow: The first is the name of the option. There 
arc six program options, each corresponding to a program variaMe. The second field is the value of the 
program variable, so the format depends on the variable. For the integMeth, minMeth, and gradMeth 
variables, the second field is a character string, left-justified in a field 16 characters wide. For the 
gradDclt and bound variables, the second field is a floating-point field. For the msgLevel variable, the 
second field is is an integer value right- justified in a field 5 characters wide. 

Examples of option records are 


option 

integ 

runge-kutta 

option 

min 

newton-raphson 

option 

gradMeth 

single-sided 

option 

gradDelt 

1 

w 

o 

o 

o 

o 

o 

o 

o 

o 

o 

t-4 

opt ion 

bound 

. 1000000000E- 

option 

msgLevel 

50 
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APPENDIX B— HELP FILES 


The online help files used by pEst are listed in the following sections as they appear in program use 
The help files detail the exact syntax of each command. The help files are alphabetically ordered. 

B.l Abort 

abort [cmd] -- exit pEst immediately 

USAGE 

abort 

DESCRIPTION 

Exits pEst without saving current program status or writing 
computed time history. 

EXAMPLES 

abort 

SEE ALSO 

quit, save, write 

KEYWORDS 

abort command, 
abort/exit pEst 

AUTHOR James Murray - NASA Dryden 
VERSION 2.1 
DATE 11/19/85 


B.2 Bound 

bound [var] -- convergence bound 
USAGE 

show bound 

set bound Convergence bound> 

DESCRIPTION 

Bound is the value used in the convergence test of the 
minimization algorithm. If the percentage change in the cost 
function between two successive iterations is less than this 
number, then convergence is declared and the program stops 
iterating . 

Convergence tests may become more complicated in future 
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versions, possibly rendering this variable obsolsts. 
Ths default is bounds. 0001 

EXAMPLES 
show bound 
set bound 0.00001 

SEE ALSO 
show, set 

KEYWORDS 

bound variable, 
convergence bound/criterion 

AUTHOR Janes Murray 

VERSION 2.1 

DATE 11/21/85 


B.3 Constant 

const [cmd] — display or set constants 
USAGE 

const <const_list> <value> 

PARAMETERS 

const.list 

List of const nanes, separated by commas or blanks. Nay 
alternatively be ’all' (for all constants). If not 
specified, it defaults to all constants. 

value 

Floating point value for all referenced constants. If this 
value is present, all referenced constants will be set to 
the specified value. If value is omitted, the constant 
values will remain unchanged. 

DESCRIPTION 

Displays or sets value of selected constants. 

EXAMPLES 

const area 200.0 
const all 

SEE ALSO 

param, flag [cmd] 


constants [var] 


KEYWORDS 

const command, 

set/change/list/ show/display constants 

AU1H0R Jama 8 Murray - NASA Drydsn 
VERSION 2.1 
DATE 3/11/86 


B.4 Constants 

constants [var] user~defined constants 
DESCRIPTION 

The following are the names and brief descriptions of the user 
routine constants currently defined in pEst. C.G positions 
are positive aft, above, and right of the reference point. 


NAME 

DESCRIPTION 

avg.qbar 

average qbar, psf 

avg.mach 

average mach 

avg.v 

average velocity, ft/sec 

avg.alpha 

average alpha, deg 

avg. theta 

average pitch attitude, deg 

avg_beta 

average beta , deg 

avg.phi 

average roll attitude , deg 

mass 

mass, slugs 

ix 

ix, slug-ft**2 

iy 

iy, slug-ft**2 

iz 

iz, slug-ft**2 

ixy 

ixy, slug-ft**2 

ixz 

ixz, slug-ft**2 

iyz 

iyz, slug-ft**2 

area 

reference area, ft**2 

span 

reference span, ft 

chord 

reference chord, ft 

xcg 

x-coordinate of eg, ft 

ycg 

y-coordinate of eg, ft 

zcg 

z-coordinate of eg, ft 

ALSO 

nst, parameters, flags, states, responses, controls, extras 


KEYWORDS 


.12 






constant names/descriptions 


AUTHOR Janas Murray 
VERSION 2.1 
DATE 11/22/85 


B.5 Controls 


controls [var] ~~ control signals 
DESCRIPTION 

Th® following are th® names and brief descriptions of the 
control signals currently defined in pEst. 


NAME 

de 

dl 

d2 

da 

dr 

d3 

d4 


DESCRIPTION 

elevacor deflection, deg 
deflection of control surface 1, deg 
deflection of control surface 2, deg 
aileron deflection, deg 
rudder deflection, deg 
deflection of control surface 3, deg 
deflection of control surface 4, deg 


SEE ALSO 

parameters, constants, flags, states, responses, extras 
KEYWORDS 

control/ input signals/names/descriptions 

AUTHOR James Murray 
VERSION 2.1 
DATE 11/22/85 


B.G Extras 

extras [var] -- extra signals 
DESCRIPTION 

The following are the names and brief descriptions of the 
extra signals currently defined in pEst. 


NAME 

qbar 

mach 


DESCRIPTION 

dynamic pressure, psf 
mach number 



alt 


altitude, ft 


SEE ALSO 

parameters, constants, flags, states, responses, controls 
KEYWORDS 

extra signals/names/descriptions 

AUTHOR James Murray 
VERSION 2.1 
DATE 11/22/85 


B.7 Flag 

flag [cmd] — display or set user-defined flags 
USAGE 

flag <flag_list> +active 

PARAMETERS 
f lag.list 

List of flag names, separated by commas or blanks. May 
alternatively be 'all* (for all flags), ’active’ (for the 
active flags), or one of several synonyms for active. If 
not specified, it defaults to all active flags. 

♦act ive (♦on , +yes , +vary , ♦enable , +t , - inact , deact , -disable , 
-no,-f ) 

Switch to set the value of the referenced flags to active or 
inactive (or on/off, etc.) If this switch is present, all 
referenced flags will be set to the specified value. If the 
switch is omitted, the flag values remain unchanged. 

DESCRIPTION 

Displays or sets status of flags. 

The flag values are boolean (true/false, etc.) There is no 
distinction between +on, +true, ^active, etc., except that 
some of these synonyms may sound more sensible than others for 
a specific flag. The interpretation of the flag values is 
solely a function of the user routines . The core pEst program 
does nothing with the flags except handle the mechanics of 
setting, displaying arid storing them. 

EXAMPLES 

flag use_avg_qBar,use_avg_mach +no 
flag all 


M 


SEE ALSO 

pax am, const [cad] 
flags [var] 

KEYWORDS 

flag command , 

set/change/list/shos/display flags 

AUTHOR James Murray - NASA Dryden 
VERSION 2.1 
DATE 3/18/86 


B.8 Flags 

flags [var] — user-defined flags 
DESCRIPTION 

The fol losing are the names and brief descriptions of the user 
routine flags currently defined in pEst. 


NAME 


DESCRIPTION 


use_avg_qbar 

use.avg.mach 

use.avg.v 

use_avg_alpha 

use_avg_ theta 

use.avg.beta 

use_avg_phi 


use average qbar in if true 

use average mach if true 

use average velocity if true 

use average alpha if true 

use average pitch attitude if true 

use average beta if true 

use average roll attitude if true 


SEE ALSO 

flag, parameters, constants, states, responses, controls, 
extras, state 


KEYWORDS 

flag names /descript ions 

AUTHOR James Murray 
VERSION 2.1 
DATE 11/22/85 


35 


•vs: 


B.9 GradDelta 


gradDelta [var] — parameter increment to use for computing 
gradients 

USAGE 

show gradDelta 

set gradDelta <value> 

DESCRIPTION 

GradDelta is the relative step size used for computing finite 
difference gradients. The step size used for each unknown 
parameter is gradDelta times the current estimate of the 
parameter . 

To avoid problems near zero, there is a fixed minimum 
parameter value in the step size computation. If a parameter 
estimate is smaller than this limit (currently .000001) , the 
step size for that parameter is gradDelta times the limit. 

The default is gradDelta* .0000001 

EXAMPLES 

show gradDelta 

set gradDelta 0.00001 

SEE ALSO 

show, set, gradMeth 
KEYWORDS 

gradDelta variable, 

delta/increment/(step size) for computing gradients, 
sensitivity matrix 

AUTHOR James Murray 
VERSION 2.1 
DATE 11/21/85 


B.10 GradMeth 

gradMeth [var] -- method of computing gradients 
USAGE 

show gradMeth 

set gradMeth <method> 


DESCRIPTION 

GradMeth is the mot hod used by the program for computing 
gradients. Two options are available. Both are finite 
difference methods. Pest has no provision for analytical 
differentiation. 

If gradHeth is set to *1* or ’s [ingle- s ided] ’ , a forward 
difference method will be used. This is the fastest method, 
requiring only n +1 integrations to compute an n-dimensional 
gradient. It sacrifices some accuracj , however. 

If gradMeth is set to »2’ or ’dfouble- sided] » , a central 
difference method will be used. This is the most accurate 
method implemented, but requires 2n integrations to compute an 
n -dimensional gradient. Since the integrations in the 
gradient computation dominate the computational time of pEst , 
the central difference method takes about twice the time of 
the forward difference method. 

For most efficient use of computer time (when it is worth the 
trouble), you might want to use forward difference computation 
most of the time, switching to central difference for the last 
few critical iterations. 

The default is gradMeth*’ single-sided’ 

EXAMPLES 

show gradMeth 
set gradMeth single 

SEE ALSO 

show, set, gradDelta 
KEYWORDS 

gradMeth variable, 
forward/central difference method, 
gradient/sensitivity computation 

AUTHOR James Murray 
VERSION 2.1 
DATE 11/21/85 


H.ll IntegMeth 

integMeth [var] -- integration method 


USAGE 


show integMeth 

set integMeth <method> 


DESCRIPTION 

IntegMeth is the method used by the program for integrating 

the equations of motion. There are currently two valid 
values. 

If integHeth is E[uler], a first order Euler method is used. 
This simple method uses only 1 state derivative evaluation and 
1 response evaluation per time point. Although fast, this 
method is quite innaccurate. It is recommended only for rough 
approximations when computer time is of critical importance. 


If integHeth is R[unge-Kutta] a 4th order Runge-Kutta method 
is used. This method requires 4 state derivative evaluations 
and 1 response evaluation per time point. Although moderately 
expensive, its accuracy is good and it is recommended for most 
applications . 

The default is integMeth* ’Runge-Kutta’ 

EXAMPLES 

show integMeth 
set integMeth runge 

SEE ALSO 
?jhov, set 

KEYWORDS 

integMeth variable, 

Euler integration method/algorithm, 

Runge-Kutta/ (Runge Kutta)/RK integration method/algorithm 

AUTHOR James Murray 
VERSION 2.1 
DATE 11/21/85 


B.12 Iterate 

iterate Land] -- perform parameter estimation iterations 
USAGE 

iterate [<niter>] [<min_meth>] 

PARAMETERS 

niter 


Maxinun number of iterations to do; if onitted or zero, the 
equations of notion will only be integrated and the cost 
function evaluated; no parameter update sill be performed. 

min.meth 

Optional specification of the nininization method to be 
used. If this parameter is included, the minMeth variable 
is set accordingly and the specified method is used. If 
this parameter is not specified the current value of the 
minMeth variable is used. The minMeth variable can also be 
set with the set command. The 3 currently valid values are 
n [ewton-raphson] g[radient] and d[avidon-fletcher-powell] . 
The initial default for min_neth is newton-raphson. See the 
minMeth helpFile for details. 

DESCRIPTION 

The iterate command initiates the parameter estimation. 
Anything beginning with ’it* eill be accepted as an 
abbreviation. 

EXAMPLES 
iterate 
iter 3 grad 

SEE ALSO 

set, show, minMeth, integMeth, gradMeth, bound, msgLevel, 
param, state, response, flag 

KEYWORDS 

itorate/iter/it command, 

set minMeth/ (minimization method) variable, 

estimate parameters, integrate equations, minimize cost 

function 

AUTHOR James Murray - NASA Dryden 

VERSION 2.1 

DATE 11/21/85 


11.13 MinMeth 

minMeth [var] -- cost function minimization method 

USAGE 

show minMeth 

set minMeth minimization method> 


DESCRIPTION 


U ** d *’ th * [ ’ r0 * r “ '° r ■ininizing 
r coat function. Thar* are currantly three valid valuta* 
nCanton-r.ph.on], gradient], and dtavidon-fletcher-powell] . 

is f ° r *>■* •itmtw ia 

actually implemented as a Gauaa-Iavton algorithm. Craaer-Rao 
bound, are available only when Mewton-Rapfson ia uaed^he 
algorithm is modified by the addition of a line search In the 
Gauss -Newton direction. With this line search, the coat 

MV * T ^ 1 “ Cr *— fr ° B «• ^-^ion to 
n ® ' “ ay not decr ®ase much (or at all), but it will 

reporterC a p 8 ; ogr ^ y bu V g 8 0litiOn ‘ * 8h ° Uld b * 

Gradient (also known as steepeat descent) is in general less 
fficient than Newton-Raphaon. The gradient algorithm is the 
mos conservative of those implemented in that ft has^ 

iliT*- ^ 18 b — *• 

involves no matrix inversions. Mote, however, that matrix 
inversion difficulties are usually a sign of more serious 
intrinsic identif lability problems and should not be lightly 
g ored. Like the pEst implementation of Gauss-Newton the 

ft. a rat"on"to h zh ‘ h °“ ld " 8, * r gi '* * incre *«« fro. 'one 

■ conwg “ c * ° f th * — i. 

it ztztvit : 0i -r Fi,td, ~ u * ig ° rui “ *■> 

. V • , 8 80 18 8ti11 somewhat experimental. The 

Ho‘^r . ' h !r r \ lly '» »• "po.erful." 

. « have not yet spent much time investigating its 

performance in the pEst environment. * * itB 

The default is minMeth-newton-raphson’ 

EXAMPLES 

show min.meth 
set min.meth newton 

SEE ALSO 

show, set, iterate 
KEYWORDS 

minMeth variable, 

»e»ton-Rzph.on/»,.ton/((;«ui,-|UM 0 „)/ (:suM method/zlzorlthn 
grzdient/Cetcepeot dncnt) nethod/elgori thn 
D.vidon-Fl.,ch.r-Po».U/dfp ..thod .Igorith.' 
co,t (unction nlnlmiz.tion/.ln/optLiz.tion ..thod/nlgorlth. 
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AUTHOR James Hurray 
VERSION 2.1 
DATE 11/21/8S 


B.14 MsgLevel 

msgLevel [var] - - Mount of screen output 
USAGE 

show msgLevel 

oat msgLevel <aa8aaga_laval> 

DESCRIPTION 

MsgLevel is tha amount of screen output displayed. The 
following describes tha output added at each value of 
msgLevel. Higher values of msgLevel also include all of the 
output for the lover levels. The numeric values assigned 
leave lots of room for future finer control. 

LEVEL DESCRIPTION 

(levels 1-9 control output that occurs only on error 
termination) 

5 Show error specifics, detailing where errors occur. 


(levels 10-49 control output that occurs once per iteration) 

10 Show cost function value each iteration. 

15 Show cost per response signal each iteration. 

20 Show estimated parameter values each iteration. 

40 Summarize each iteration. 

(levels 50-99 mostly control output that can be many times 

per iteration) 

SO Summarize each line search. 

52 Summarize each stage of line search. 

55 Show all integration failures, recoverable or not. At 
lower message levels, some integration failures are 
silently recovered, for instance by trying smaller step 
sizes . 
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57 Show correlation matrix every newton-raphson iteration. 

60 Detail each cost evaluation in line searches. This 
tracks the convergence of the line search algorithm. 
Mostly used in development of line search algorithms. 

62 Detail integration failures. Show details of why 
integration was terminated. 

65 Show second gradient matrices every newton-raphson 
iteration. Show approximate inverse second gradient 
matrices every dfp iteration. 

70 Show all gradient vectors. 

75 Show cost function value whenever it is evaluated for 
any purpose. 

80 Show the cost per response signal whenever the cost is 
evaluated . 

(levels 100-? control output that occurs every time point of 

integration) 

(Note that the volume of output can be extremely large) 

100 Show the state vector every time point. (Unimplemented) 
The default is msgLevel s 50 

EXAMPLES 

show msgLevel 
set msgLevel 60 

SEE ALSO 
set , show 

KEYWORDS 

msgLevel variable, 
msg, 

level/amount of output messages 

AUTHOR James Murray 

VERSION 2.1 

DATE 11/22/85 
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B.15 Parameter 


paraa [cad] -- display or set parameter values or status 

USAGE 

par am <param_list> <value> 

♦active ♦reset 

♦predicted ♦crBound ♦delta 

PARAMETERS 

param.list 

List of parameter names, separated by commas or blanks. May 
alternatively be ’all' (for all parameters), ’active’ (for 
the active parameters), or one of several synonyms for 
active. If not specified, it defaults to all active 
parameters. In this context, an ’active’ parameter is one 
being allowed to vary in the estimation (varying is one of 
the acceptable synonyms) . 

value 

Floating point value for all referenced parameters. If this 
value is present, all referenced parameters vill be set to 
the specified value. If value is omitted, the parameter 
values vill remain unchanged. 

♦active (+on, ♦yes, ♦vary,+enable,+t,-inact,-deact, -disable, 
-no,-f) 

Switch to activate or deactivate estimation of all 
referenced parameters. If ^active is set, the referenced 
parameters will be estimated. If -active is set, they will 
not be estimated. If neither is set, the status of the 
parameters will remain unchanged. 

♦reset 

If ♦reset is specified, all parameters in the param.lipt 
will be reset to their predicted values. Mote that it does 
not make much sense to specify both a value and ♦reset; if 
you do so, the Preset specification will override without 
comment. The switch may be abbreviated to +r. The default 
is -reset. 

♦predicted 

If the ♦predicted switch is specified, the predicted values 
of the parameters will be included in the display resulting 
from this command. The switch may be abbreviated to +pr. 

The default is -predicted. 


♦crBound 


If the ♦crBound switch is specified, the Cramer-Rao bounds 
of the parameters sill be included in the display resulting 
from this command. Anything beginning with +cr or +b is 
accepted as a synonym. The default is -crBound. 

♦delta 

If the ♦delta switch is specified, the the parameter value 
changes in the previous iteration will be included in the 
display. Anything begining with +d will be accepted as a 
synonym (unless the second letter is i, which could cause 
confusion with +disable). 

DESCRIPTION 

Displays or sets value and estimation status of selected 

parameters. Any command starting with 'par' will be accepted 
as a synonym. 

EXAMPLES 

param cma.cmde.cmq +fix 
param cnp -0.25 
parm vary 
par all 

SEE ALSO 

const, flag [cmd] 
parameters [var] 

KEYWORDS 

param command, 

set/change/list/show/display parameter /param/parm values/ 
status 

AUTHOR James Murray - NASA Dryden 

VERSION 2.1 

DATE 3/18/86 


B.1G Parameters 

parameters [var] -- estimated parameters 
DESCRIPTION 

The following aie the names and brief descriptions of the 
estimated parameters currently defined in pEst. Instrument 

positions are positive aft. above, and right of the reference 

point . 


U 


DESCRIPTION 


NAME 

cNormO 

cNorma 

cNormq 

cNormde 

cNormdl 

cNormd2 

cNormd3 

cmO 

cma 

cmq 

cmd© 

cmdl 

cmd 2 

cmd3 

caO 

caa 

caq 

cade 

cadi 

cad2 

cad3 

vO 

alphaO 

qO 

thetaO 

vBias 

alphaBias 

qBias 

thetaBias 

anBias 

axBias 

qdotBias 

ka 

xa 

ya 

za 

xan 

yan 

zan 


normal force aerodynamic bias 
normal fore* due to alpha, par deg 
normal force due to pitch rate, per rad 
normal force due to elevator deflection, per 
deg 

normal force due to dl deflection, per deg 
normal force due to d2 deflection, per deg 
normal force due to d3 deflection, per deg 
pitching moment aerodynamic bias 
pitching moment due to alpha, per deg 
pitching moment due to pitch rate, per rad 
pitching moment due to elevator deflection, 
per deg 

pitching moment due to dl deflection, per deg 
pitching moment due to d2 deflection, per deg 
pitching moment due to d3 deflection, per deg 
axial force aerodynamic bias 
axial force due to alpha, per deg 

force due to pitch rate, per rad 
axial force due to elevator deflection, per 
deg 

axial force due to dl deflection, per deg 
axial force due to d2 defection, per deg 
force due to d3 deflection, per deg 
increment to velocity initial condition, 
ft/sec 

increment to alpha initial condition, deg 
increment to pitch rate initial condition, 
deg/sec 

increment to pitch attitude initial condition 
deg 

measurement bias on velocity, ft/sec 
measurement bias on alpha, deg 
measurement bias on pitch rate, deg/sec 
measurement bias on pitch attitude, deg 
measurement bias on normal acceleration, g 
measurement bias on axial acceleration, g 
measurement bias on pitch acceleration, 
deg/sec**2 

upvash factor for alpha sensor 
x-coordinate of alpha sensor, ft 
y-coordinate of alpha sensor, ft 
z-coordinate of alpha sensor, ft 
x-coordinate of normal accelerometer, ft 
y-coordinate of normal accelerometer, ft 
z-coordinate of normal accelerometer, ft 




xax 

yax 

zax 

xv 

yv 

zv 

cyO 

cyb 

cyp 

cyr 

cyda 

cydr 

cyd3 

cyd4 

clO 

clb 

clp 

clr 

clda 

cldr 

cld3 

cld4 

cnO 

cnb 

cnp 

cnr 

cnda 

cndr 

cnd3 

cnd4 

betaO 

pO 

rO 

phiO 

betaBias 

pBias 

rBias 

phiBias 

ayBias 

pdotBias 

rdotBias 


x-coordinate of axial accelerometer, ft 
z-coordinate of axial accelerometer , ft 
z-coordinate of axial accelerometer, ft 
x-coordinate of velocity sensor , ft 
y-coordinate of velocity sensor, ft 
z-coordinate of velocity sensor, ft 
side force aerodynamic bias 
side force due to beta, per deg 
side force due to roll rate, per rad 
side force due to yaw rate, per rad 
side force due to aileron deflection, per deg 
side force due to rudder deflection, per deg 
side force due to d3 deflection, per deg 
side force due to d4 deflection, per deg 
rolling moment aerodynamic bias 
rolling moment due to beta, per deg 
rolling moment due to roll rate, per rad 
rolling moment due to yaw rate, per rad 
rolling moment due to aileron deflection, per 
deg 

r °Hing moment due to rudder deflection, per 
deg 

rolling moment due to d3 deflection, per deg 
rolling moment due to d4 deflection, per deg 
yawing moment aerodynamic bias 
yawing moment due to beta , per deg 
yawing moment due to roll rate, per rad 
yawing moment due to yaw rate, per rad 

yawing moment due to aileron deflection, per 
deg 

yawing moment due to rudder deflection, per 
deg 

yaving moment due to d3 deflection, per deg 
yawing moment due to d4 deflection, per deg 
increment to beta initial condition, deg 
increment to roll rate initial condition, 
deg/sec 

increment to yaw rate initial condition, 
deg/sec 

increment to bank angle initial condition, deg 
measurement bias on beta, deg 
measurement bias on roll rate, deg/sec 
measurement bias on yaw rate, deg/sec 
measurement bias on bank angle, deg 
measurement bias on lateral accelerometer, g 
measurement bias on roll acceleration, 
deg/sec**2 

measurement bias on yaw acceleration. 


•Hi 


* 
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kb 

deg/sec**2 

sidevash factor for beta sensor 

xb 

x- coordinate of beta sensor, ft 

yb 

y-coordinate of beta sensor, ft 

zb 

z-coordinate of beta sensor, ft 

xay 

x-coordinate of lateral accelerometer 

yay 

y-coordinate of lateral accelerometer 

zay 

z-coordinate of lateral accelerometer 

gAlpha 

feedback gain for alpha state 


feedback gain for q state 

gBeta 

feedback gain for beta state 

gP 

feedback gain for p state 

g R 

feedback gain for r state 


SEE ALSO 

par an, constants, flags, statss, responses, controls, extras 
KEYWORDS 

parameter/param names/descriptions 

AUTHOR James Murray 
VERSION 2.1 
DATE 3/10/86 


B.17 pEst 

pEst [cmd] -- parameter estimation program. 

USAGE 

C/user /murray /consnands/] pEst 
DESCRIPTION 

This program does parmater estimation. The program is 
designed for interactive operation. 

The program resides in the directory /us er /murray /pEst , with 
aliases .in /user/murray/commands. The usePest command 
facilitates access to pEst . 

There is a full internal help facility which covers the 
commands within pEst. 

EXAMPLES 

/ user/murray/commands/usePest 
pEst 

read measured 
restore current 
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iter 3 newton 

plot 

quit 

CAVEATS 

Current dimensions limit the program to 200 parameters and 
2000 time points. The limits are all checked, so that 
accidentally exceeding them sill not cause the program to 
crash. 

ERROR HANDLING 

The program attempts to recover from all errors. Such mundane 
errors as exceeding dimension limits, or giving a non-existent 
file name or signal name are all caught. The program should 
not crash, regardless of what junk you feed it for commands. 
Infinite or NaN quantities in the data may crash it. If you 
succeed in crashing the program in any other way, please let 
me know. 

SEE ALSO 

bindPEst, thPlot, usePest, internal help 

FILES 

current current program status. 

measured time history data file. 

computed time history file of estimated variables. 

pEst_thplot .temp scratch file for communicating with thplot 

IMPLEMENTATION 
Fortran program. 

The time history data file interface routines are used to read 
the data files. See the help topic f ilelnterf ace for 
discussions of the file interface subroutines. You must write 
customized versions of these routines to use pEst on data 
files not supported by the standard ones. 

KEYWORDS 

program pEst, 
parameter estimation, 
mmle 

AUTHOR James Murray - NASA Dryden 
VERSION 2.1 
DATE 10/28/85 
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B.18 Plot 


plot [cffld] plot fit* of currently active responses, control, 
and extras 

USAGE 

plot C+resCiduals]] [+ u I *c Controls] I ♦iCnputs]] [+ex[tras]] 

PARAMETERS 

♦res 

If this flag is turned on, the fit residuals are also 
plotted. Default is -res. 

♦u 

If this flag is turned on, the contol signals are also 
plotted. Default is -u. 

+ex 

If this flag is turned on, the extra signals are also 
plotted. Default is -ex. 

DESCRIPTION 

Enters thPlot program and automatically plots time history 
fits of all currently active response variables. Fit 
residuals and other signals are selectively plotted. Upon 

completion of the time history plots, control is returned to 
pEst . 

This command differs from the thplot command in that the 
thplot command expects you to interactively enter any commands 
to the thplot program. That allovs more flexibility in what 
is plotted, but requires more input from the user. The plot 
command, in contrast neither expects nor allows interactive 
input during the plotting. 

EXAMPLES 

plot 

plot +res 
plot -*u ^extra 

CAVEATS 

As currently implemented, the plot command creates a file 
called pEsc. thplot .temp, used for input to the thplot program. 
Any pre-existing user file of that name will be overwritten 
(and later deleted). 

ERROR HANDLING 

Automatically uses the write command to create the computed 


t 
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tine history file. Therefore, all of the error handling 
discussed under that command applies. If the write command 
fails to write a computed time history file, the plot command 
will terminate and return to the pEst command line without 
doing any plots. If a bad flag is specified, the plot command 
will terminate and return to the pEst command line without 
doing any plots. 

SEE ALSO 

thplot, write 

KEYWORDS 

plot command , 

plot response fits, residuals, controls/ inputs , and extras 

AUTHOR James Murray - NASA Dryden 

VERSION 2.1 

DATE 11/19/85 


B.19 Quit 

quit [and] — save current status and exit pEst 

USAGE 

quit 

DESCRIPTION 

Writes current program status and computed time history files, 
and exits pEst. The effect is the same as doing a save 
command with no parameter, followed by an abort command. 

EXAMPLES 

quit 

SEE ALSO 

abort, write, save 

KEYWORDS 

quit command, 

save current status and quit/exit pEst 

AUTHOR James Murray - NASA Dryden 
VERSION 2.1 
DATE 11/19/85 
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B.20 Read 

read [and] — read measured time history iile 
USAGE 

read [fileName] 

PARAMETERS 

filename 

filename (or pathname) of the file to be read. If not 
specified, it defaults to the same name as last specified on 
a read command. If no filename has been specified by any 
read command in the current run, the name defaults to 
’measured’ . 

DESCRIPTION 

Reads a measured time history data file. A measured time 
history file must be read before pEst can do anything of 
consequence. On initial startup, the program attempts to read 
a time history file named ’measured’ . If this attempt fails 
for any reason, you will be unable to do much until you do a 
successful read command. The read command can also be used, 
after completing analysis of one maneuver , to switch to a new 
data set for analysis. 

A few program variables are automatically reset whenever a 
read command is successfully executed. These are variables 
that are closely associated with the individual maneuver and 
unlikely to be meaningfully carried over from one case to 
another . 

Primarily, all maneuver and window times are redefined. The 
maneuver times are automatically determined by looking for 
time dropouts of greater than 1 second (exact number subject 
to future change). The window times are set to equal the 
maneuver times by default. 

The avg_ constants are also reset to the average of the 
measured data for the maneuver. It is possible to manually 
override these values if you want to (for instance if the 
measured values are wrong) , but it is unlikely that you want 
such overrides to automatically carry over from one maneuver 
to another (in fact it could cause great confusion if you were 
not aware that such carryover was occuring). Therefore, 
average values are reset each time a maneuver is read. 

EXAMPLES 

read 
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read measured. casel3 


ERROR HANDLING 

The program gives error messages and returns to the command 
line on encountering any error. Previously read data may be 
overwritten by a failed attempt to read new data; in this 
case, the program will clear out all of the corrupted data and 
you will have to execute a new successful read couand before 
doing much else constructive . 


All commands that use time history data will recognize and 
appropriately handle the case where no valid time history data 
is available. (Appropriate handling usually means that the 
command will give an error message and refuse to execute) . 

SEE ALSO 

write, save, restore 

KEYWORDS 

read command, 

read/load time history data file 

AUTHOR James Murray 
VERSION 2.1 
DATE 11/19/85 


11.21 Response 

response/output [cmd] -- display or set response equation status 
USAGE 

response/output <resp_list> ^active weight«<weighting> 
resp.list 

List of response variable names, separated by commas or 
blanks. Hay alternatively be *all’ (for all responses), 
’active’ (for the active responses), or one of several 
synonyms for active. If not specified, it defaults to all 
active responses. 

♦active (♦on.+yes,tvary,4enable,+t.-inact,-deact, -disable, 

-no,-f ) 

Switch to activate or deactivate computation of the 
referenced responses. If ^active is set, the referenced 
response equation(s) will be activated. If -active is set, 
they will be deactivated. If neither is set, the status of 
the response equations will remain unchanged. 


r »2 
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weighting 

Floating point value for response error weighting. If this 
parameter is specified, the cost function weightings of all 
the referenced outputs will be set to the specified value. 

If it is not specified, the weightings will be left 
unchanged. The keyword may be abbreviated to anything 
starting with w. 

DESCRIPTION 

Displays or sets status of response variable(s). Anything 
beginning with ’resp’ or ’out’ will be accepted for this 
command . 

Note that the weighting is irrelevant for responses that are 
inactive. An inactive response makes no contribution to the 
cost function, making its effective weighting zero. The 
differences between setting the weighting to 0 and 
deactivating the response are two: First, the responses are 
still computed (just not used in the cost function) if the 
weighting is zero for an active response; if the response is 
inactive, it is never computed. Obviously, if you need the 
response computation for some reason other than its use in the 
cost function (perhaps you want a plot of it) , you must have 
the response active. This may involve substantial computation 
time, so if you don’t need the response, it is pref err able to 
deactivate it. The second difference is just one of 
convenience: if you deactivate a response, the program still 
remembers what its weighting was in case you later want to 
re-activate it with the same weighting as before. 

EXAMPLES 
response 
resp alpha +off 
output beta ^activate w*100.0 

SEE ALSO 
state [cmd] 
responses [var] 

KEYWORDS 

response command, 

set/list/show/display response/output equation/variable 
status/weighting 


AUTHOR James Hurray - NASA Dryden 
VERSION 2.1 
DATE 3/18/86 


B.22 Responses 

responses [var] -- response signals 
DESCRIPTION 

The following are the naaes and brief descriptions of the 


response 

signals currently defined in pEst. 

NAmE 

DESCRIPTION 

V 

velocity, ft/sec 

alpha 

angle of attack, deg 

q 

pitch rate, deg/sec 

theta 

pitch attitude, deg 

an 

normal acceleration, g’s 

ax 

longitudinal acceleration, g’s 

qdot 

pitch acceleration, deg/sec**2 

beta 

angle of sideslip, deg 

P 

roll rate, deg/ sec 

r 

yaw rate, deg/sec 

phi 

bank angle, deg 


lateral acceleration, g’s 

pdot 

roll acceleration, deg/sec**2 

rdot 

yaw acceleration, deg/sec**2 

SEE ALSO 



response, parameters, constants, flags, states, controls, 
extras 

KEYWORDS 

response/ output signals/names/descriptions 

AUTHOR James Murray 
VERSION 2.1 
DATE 11/22/85 


11.23 Restore 

restore [cmd] -- restore status from file 

USAGE 

restore [fileName] 

PARAMETERS 
f ilename 

filename (or pathname) of the file to be read. If not 
specified, it defaults t;> the same name as last specified on 
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a restore command. If no fileName has been specified by any 
restore command in the current run, the name defaults to 
’current* . 

DESCRIPTION 

Reads program status from the specified file. The program 
status includes the values of all parameters and program 
options. The file can be from either of two sources. The 
file may have been created by an independent program, usually 
by looking up predicted derivative values from a data table; 
this is often the case when starting a new case. 

Alternatively, the file can have been created from a save 
command; in this case, the program will be restored to the 
same status as when the save command was executed. 

A file created while analyzing one case can be restored while 
analyzing a different set of time history data. This allows, 
for instance, a convenient way of starting the estimation from 
converged values obtained from a case at a similar flight 
condition. Some values, such as the parameter estimates, can 
be meaningfully transfered from one case to another in this 
manner. However, other values on the current status file are 
unlikely to be useful (or legal) when applied to a different 
time history data file. 

In particular, the window definitions are closely tied to 
specific maneuvers. Window times for one case are likely to 
be outside of the maneuver times for a different case. 

The avg_ constants are also closely tied to the specific 
maneuvers ; they usually equal the average measured values for 
the maneuver. It is possible to manually override these 
values if you want to (for instance if the measured values are 
wrong) , but it is unlikely that you want such overrides to 
automatically carry over from one maneuver to another (in fact 
it could cause great confusion if you were not aware that such 
carryover was occuring) . 

To avoid these undesired carryovers, while allowing carryover 
of other pertinent parameters, the restore command checks for 
consistency between the status file and the time history data 
being analyzed. The status file includes a record of the 
maneuver times and average measured signal values to 
facilitate this comparison. If the data on the status file is 
consistent with the time history data being analyzed, then all 
parameters from the status file are accepted, including window 
times and avg_ constants. If the data on the status file is 
inconsistent with the maneuver being analyzed, then the window 




tines and avg_ constants on the status file are rejected; 
instead, the windows are set equal to the available Maneuvers 
and the avg_ constants are set equal to the average measured 

values . 

For similar reasons, the window times and avg_ constant values 
are reset every time a read command is executed. 

On initial startup, if the automatic read command works, the 
program automatically attempts a restore command for the file 
name ’current’. 

Unlike the read command, which must succeed before you can 
analyze any data, the restore command is technically optional. 
It is possible to analyze a case without ever reading a status 
file. All of the it^ms set by the restore command have 
default values. In some cases where no meaningful default is 
possible, the default value is an illegal one (for instance 
weight defaults to 0) . Such values must be reset either by a 
restore command or by manual entry before estimation can 
procede. The program will give an error message if you 
attempt estimation without setting these variables to legal 
values . 

EXAMPLES 

restore 

restore current . casel3 
CAVEATS 

Because of the interaction between the read and restore 
commands, you must be cautious of the order in which you 
execute them. If you want to do a read command and a restore 
command to completely restore a previously saved status, the 
read command must precede the restore command. 

The data on the saved file is in ascii form, accurate to 
’only’ about 10 digits. For ’reasonable’ situations, this 
should be far more accuracy than you will need. For highly 
unstable systems, it is possible for the rounding to 10 digits 
to make the restored state noticably different than that 
originally saved. The validity of the estimates for any case 
in which this effect is noticable is in serious doubt anyway. 

SEE ALSO 
save, read 

KEYWORDS 

restore command, 
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restore/read/load currant status 


AUTHOR James Hurray 
VERSION 2.1 
DATE 11/19/85 


B.24 Save 

save [cmd] -- save current status to a file. 

USAGE 

save [fileName] 

PARAMETERS 

filename 

filename (or pathname) of the file to be written. If not 
specified, it defaults to the same name as last specified on 
a RESTORE command. If no fileName has been specified by any 
restore command in the current run, the name defaults to 
’ current * . 

DESCRIPTION 

Saves current program status to a file. The program status 
includes the values of all parameters and program options. 

The resulting file can be used for several purposes: The 
program can later be restored to the saved point by using the 
restore command. The file can be input to plotting programs 
or other programs that use the parameter values estimated by 
pEst. The file can also be printed (it is ascii) as a record 
of the results. 

EXAMPLES 

save 

save current. case23. 19Nov85 
CAVEATS 

Note that the default fileName is taken from the last RESTORE 
command, rather than the last SAVE command. The logic behind 
this is that a save command with no arguments is taken as a 
request to save the current status in place of the file that 
you started with; this is the most usual mode of operation. A 
save specifying a different fileName is assumed to be asking 
to save a specific status to a specific file for special 
purposes. All subsequent work reverts to the original file. 

DUGS 

By normal Fortran defaults, the file is created with the 


♦fortranCCTL flag sat, in spita of tha fact that tha file doas 
not have fortran carriage control characters in column 1. In 
order to correctly print the file, you must explicitly specify 
-fortranCCTL in tha print command (or do a modifyFile to 
correctly sat tha flag Maintained with tha file) . 

SEE ALSO 

restore, write 

KEYWORDS 

save command, 
write/save current status 

AUTHOR Jamas Hurray 

VERSION 2.1 

DATE 11/19/85 


n.25 Set 

set [cmd] — set program variables. 

USAGE 

set <variable-name> <value> 

PARAMETERS 

variable-name 

Name of tha variable to be set. 
value 

Value to be used for the variable. Tha syntax and legal 
values may vary for different variables. In particular, 
some variables have several components and the value syntax 
must specify which component is being set. 

DESCRIPTION 

The set command sets the values of program variables. For a 
list of the available variables, do "help variables". For 
details about a specific variable, do help on the variable 
name. 

EXAMPLES 

set msgLevel 80 
set integMeth euler 
set gradDelta 0.00001 
set vindow time 0.6 9.5 

SEE ALSO 


58 


show 


KEYWORDS 

set command, 

set/change program variable values 

AUTHOR James Hurray - NASA Dryden 
VERSION 2.1 
DATE 11/20/85 


B. 26 Show 

show/list/display [cmd] — show values of program variables 
USAGE 

shov/list/display <variable-name> 

DESCRIPTION 

Shows the value of the specified program variable. For a list 
of the available variables, do "help variables". For details 
about a specific variable, do help on the variable name. 

List and display are accepted as synonyms. Sh and disp are 
acceptable abbreviations. 

EXAMPLES 
list integ 
show bound 
display window 
show statistics 

SEE ALSO 
set 

KEYWORDS 

show/list/display command, 
show/list/display program variable values 

AUTHOR James Hurray - NASA Dryden 
VERSION 2.1 
DATE 11/20/85 


I). 27 Slat* 

state [cmd] -- display or set state equation status 
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USAGE 

stats <state_list> +active limit»<limit> 

PARAMETERS 

state.list 

List of stats names, ssparatsd by commas or blanks. May 
altsrr.ativsly bs 'all* (for all statss) , ’active’ (for ths 
activs statss), or ons of ssvsral synonyms for activs. If 
not spscifisd, it dsfaults to all activs statss. 

♦active ( +on , ♦yss , +vary , ^enable , +t , - inact , -dsact , -disabls , 
-no,-f) 

Switch to activats or dsactivats intsgration of ths 
rsfsrsncsd statss. If ♦active is sst, ths refsrsncsd stats 
squation(s) will bs intsgratsd. If -activs is sst, they 
will not bs integrated. If neither is set, the status of 
the stats equations will remain unchanged. 

limit 

The limit on ths absolute value of the referenced statss 
allowed during integration. Any intsgration that exceeds 
this limit will bs abandoned. If this parameter is not 
specified, the previous limits remain unchanged. 

DESCRIPTION 

Displays or sets status of state equation(s). 

The system equations can change fairly significantly depending 
on which states are active and which ones are inactive. If 
the state is inactive, no computed value for that state 
variable is defined; thus response equations or other state 
equations that use the deleted state variable must be revised. 
The revision done is to substitute some other quantity for the 
unavailable state variable. 

By default, if a state is deactivated, the corresponding 
measured quantity is substituted. Any output bias on the 
measured quantity is subtracted before substituting it. Also 
the measured flow angles are corrected to the center of 
gravity before substituting them. 

The variables with names like use.avg.alpha can force 
alternate substitutions. If use.avg.alpha is true, then the 
constant avg.alpha is substituted wherever the computed alpha 
state would otherwise have been used. In this case, the 
substitution occurs even if the alpha state is active. 
Corresponding comments apply to the other states. 


EXAMPLES 

state thota +on limit- 1000 

state all 

state phi -active 

SEE ALSO 

response [end] 
states [var] 

KEYWORDS 

state command, 

set/change/list/shoe/display state equation/variable status, 
add/delete state equations 

AUTHOR James Murray - NASA Dryden 
VERSION 2.1 
DATE 3/17/86 


D.28 States 

states [var] — state variables 
DESCRIPTION 

The folloeing are the names and brief descriptions of the 
state variables currently defined in pEst. 


MAKE 

DESCRIPTION 

V 

velocity, ft/sec 

alpha 

angle of attack, deg 

q 

pitch rate, deg/sec 

theta 

pitch attitude, deg 

beta 

angle of sideslip, deg 

P 

roll rate, deg/sec 

r 

yaw rate, deg/sec 

phi 

bank angle, deg 


SEE ALSO 

state, parameters, constants, flags, responses, controls, extras 
KEYWORDS 

state variables/naaes/descriptions 

AUTHOR Janes Murray 
VERS I UN 2.1 
DATE 11/22/85 


D.29 Statistics 


stats [var] — time history statistics 
USAGE 

show statistics] 

DESCRIPTION 

The 'shoe statistics' command lists the average value of every 
measured response, input, and extra signal for the time 
history data currently loaded. 

Statistics is not really a 'variable' as it can not be set; it 
can only be displayed. 

EXAMPLES 
show stats 
show statistics 

SEE ALSO 
show 

KEYWORDS 

statistics variable, 

maneuver/signal statistics/stats/averages 

AUTHOR James Murray 
VERSION 2.1 
DATE 11/21/85 
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thplot [cad] -- execute thPlot program 

USAGE 

thplot 

DESCRIPTION 

Writes time histories of computed responses, then enters 
thPlot program. Upon termination of the thPlot program, 
control is returned to pEst. The computed responses are 
written by an automatic write command with no arguments; 
therefore the file name is vhatever such a write command would 
use (usually the name ’computed*) 

This command differs from the plot command in that the thplot 
command expects you to interactively enter any commands to the 
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thplot program. This allows mors flsxibility in what is 
plotted, but requires more input from the user. The plot 
command, in contrast neither expects nor allows interactive 
input during the plotting. 

EXAMPLES 

thplot 

ERROR HANDLING 

Automatically uses the write command to create the computed 
time history file. Therefore, all of the error handling 
discussed under that command applies. If the write command 
fails to write a computed time history file, the plot command 
will terminate and return to the pEst command line without 
doing any plots. 

SEE ALSO 
plot, write 

KEYWORDS 

thplot command, 
time history plots 

AUTHOR James Murray - NASA Dry den 

VERSION 2.1 

DATE 11/19/8S 


B.31 Title 

title [var] — title for current file and plots. 

USAGE 

show title 

set title ’<title>’ 

DESCRIPTION 

Title is the title used on the current file and on the top of 
each plot page. The length of the title is limited to 40 
characters; longer titles are truncated without comment. If 
the title contains imbedded blanks (not uncommon) , the title 
must be enclosed in quotes; otherwise the title will be 
truncated to the first word. 

The default is title ■ » * . 

EXAMPLES 
show title 


s«t title ’ x29 flight 20 case 10* 


SEE ALSO 
8 how , sat 

KEYWORDS 

titla variable, 
plot titla 

AUTHOR James Murray 
VERSION 2.1 
DATE 2/19/86 


B.32 Version 

version [topic] — changes in pEst with version 2.1 
DESCRIPTION 

pEst version 2.1 is now released. The following describes the 
changes between version 2.1 and the previous release, version 
2.0. Access to the new version is automatic when you use pEst 
normally. The previous version will be retained for an 
interim period as a backup. To access the previous version, 
use the command pEst2.0. 

COMPATABILITY 

Version 2.1 has several small differences from version 2.0. 

The exact syntax of a few of the commands has changed. This 
require getting used to the revised syntax and changing 
any command files; the changes are all small. 

If you have a customized set of user routines, you must make a 
minor change to the argument list of defineNames; this change 
is small and mostly involves deleting a few lines from the 
routine. 

There are several additions to the current file, but version 
any current files from version 1.6 or later will be accepted 
without error. Any program that reads current files must be 
able to accept the new fields. 

USER ROUTINES 

All the arguments have been removed from subroutine 
defineNames. They were not being used for anything 
constructive anyway and they caused excessive clutter in the 
main program. 
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The see antics of the interface to subroutine computeZ have 
changed. Wa nos alios computeZ to change its xc (computed 
state) argument. Any such change constitutes the correction 
step of a predictor-corrector algorithm such as a Kalman 
filter. This change is solely one of establishing an expanded 
convention. It requires no code changes in old code, but just 
liberalizes shat nes code is allosed to do. 

USAGE 

The par am command has been expanded to allow optional display 
of or reset to the predicted values. See the param helpFile 
for details. 

The title variable has been added to allow the title used on 
the current file to be modified. The sane title is now also 
put on time history plots. See the title helpFile for 
details. 

The plot command uses thPlot's +expand option to expand scales 
to the full screen width. 

An extension to the state command allows the state variable 
limits to be displayed and changed. These limits are also now 
saved on the current file. The state variable limits are what 
stops the integration when it goes unstable. 

There are S nes parameters (gAlpha,gQ,gBeta,gP,gR) in the 
default user routines. These parameters are gains in the 
filter error formulation. Among other things, they can be 
used to stabilize the estimation for unstable systems. Their 
use at this time is still experimental. 

DOCUMENT AT I OH 

The helpfiles mentioned in the following "SEE ALSO" section 
have been modified to reflect the changes. The helpFile 
format for showing command syntax has been redone to be more 
consistent with Elxsi helpFiles (also a little simpler) . This 
file is included under "help version". 

SEE ALSO 

param, const .flag, state, response, 

KEYWORDS 

pEst version 2.1 changes 

.VJTKOR James Murray 
VERSION 2.1 
DATE 3/17/86 


B-33 Window 


window [var] — maneuver window descriptor 
USAGE 

show window 

set window [<window number>] [man[*] <maneuver number>] 
[time[=]<start> <end>] 

DESCRIPTION 

A window is a time interval which is wholly contained within a 
maneuver time interval. A maneuver may contain multiple 
windows . 

When the equations of motion are evaluated, the computed time 
history in generated only for those time points contained 
within the window(s) . Outside of the window(s), the average 
value for the output variables is used. Also, the cost 
function for the estimation process depends only on those time 
points within the window(s). 

If, when a window is defined, the window number is not 
specified, all existing windows will be deleted and the new 
window will be window 1. The window number, if specified, 
cannot exceed the current window count by more than 1. If the 
maneuver number is rot specified, the default maneuver number 
1 will be used. The window start and end times must be 
specified and are given relative to the starting time of the 
maneuver. 

The default is that the windows exactly equal the available 
maneuvers . 

EXAMPLES 

show windows 

set wind 1 time*0.5-9.5 

set wind man*3 time 1.25 5.55 

SEE ALSO 
show, set 

KEYWORDS 

window variable, 
multiple maneuvers, 
start/stop times, 
time segments/intervals/windows 

AUTHOR James Murray 
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VERSION 2.1 
DATE 11/21/85 


B.34 Write 

write [end] — write computed time history file 
USAGE 

write [fileName] 

PARAMETERS 

filename 

filename (or pathname) of the file to be written. If not 
specified, it defaults to the same name as last specified on 
a write command. If no fileName has been specified by any 
write command in the current run, the name defaults to 
’computed* . 

DESCRIPTION 

Writes a computed time history data file. This file is most 
often used for plotting, but can have numerous applications. 

It could be input to programs to analyze the residual 
statistics. It also allows pEst to be used as a simple batch 
simulator for other purposes. 

The signals written to the file are the calculated states 
(-s-hat suffixes), the calculated responses (-hat suffixes), 
and the residuals (-res suffixes) . The names of all signals 
are derived by appending the indicated suffixes to the basic 
state or response variable name. The residual is defined to 
be the measured response minus the calculated response. 

EXAMPLES 

write 

write computed. casel2 
ERROR HANDLING 

If the output file can not be opened or if integration of the 
equations is disallowed for any reason, the write command will 
give an error message and return to the command line. 

If the integration fails part may through (usually because it 
goes unstable) , the file is written normally up to the point 
where the integration failed. The remainder of the file is 
written with average measured values substituting for the 
unavailable calculated responses. This is done to accomodate 
the plotting program, which wants the same time points in both 
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the measured and calculated data files in order to plot a fit. 
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SEE ALSO 

read, plot, thplot 
KEYWORDS 

vrite/save command, 

write/save computed time history data 

AUTHOR James Murray ! 

VERSION 2.1 
DATE 11/19/85 
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